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ABSTRACT 


This thesis proposes a framework by which optimal system renewal decisions 
may be made in a consistent and timely manner within the military, under realistic 
conditions of changing environments. A unique network representation of the system 
renewal process was used to develop a prototype version of the analytical core model. 
Its plausibility and usefulness was demonstrated by a series of case studies. The case 
studies also show how pertinent staff forecasts can be organised and integrated to 
provide decision makers with a broad, consistent, long-term perspective of the issues 
relevant to system renewal planning. They are presented with a graphic picture of 
the entire solution space as structured by the scenarios considered. Various solutions 
are suggested for each scenario and their robustness may be tested by thoroughly 
exercising the model for a wide range of scenarios. Prediction of what and when 
renewals are likely to be and estimation of the associated cost and effectiveness 
allow anticipatory long-term plans to be formulated. The thesis also suggests how 
bulk versus phased procurement decisions and force level and mix issues could be 


analysed using the model. 
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I. INTRODUCTION 


A. INTRODUCTION 

Military forces continually undergo change as obsolete or uneconomical assets 
are retired, new weapons are introduced and existing ones modified or upgraded. 
This process of system renewal is affected by a myriad of economic, technological 
and politico-military considerations. Changes in enemy capabilities or doctrine, for 
instance, may render certain systems less effective or even obsolete. New technology 
may offer opportunities to improve present capabilities. Most of the available open 
literature, analyses and models for system renewal are either one-dimensional, too 
general, or too commercially oriented to be useful to military planners. There is a 
need to synthesize a coherent, balanced view relevant to military planners to provide 


comprehensive policy guidance on upgrading, replacement and related matters. 


B. PROBLEM 

The system renewal planning problem is encountered prior to the last milestone 
of the defence acquisition cycle. It is basically one of being able to evaluate and 
decide in a timely and rational manner, under conditions of changing technological, 


economic and politico-military environments, both the following: 


1. When will a defence asset, in its present form, become obsolete or no longer 


cost-effective to retain (relative to existing or emerging alternatives)? 


2. What should then be done? What are the alternatives and which one should 


be chosen? 


C. DEFINITION 

System renewal encompasses all changes made to maintain or improve the 
cost-effectiveness of the present defence system. It includes the acquisition of new 
capabilities, retirement from service of obsolete ones, relegation to reserves of sur- 
plus capabilities and conversion of systems to new roles. This view streamlines the 
treatment of changes to the status quo and involves a higher level of system repre- 
sentation. It allows us to consider complex changes, such as when a new support 
ship displaces a couple of smaller ones, one of which is preserved for the boneyard 


and the other converted to a hospital ship. 


D. PURPOSE 

The purpose of this thesis is to develop a rational and consistent decision 
model for system renewal planning that is relevant and useful to military planners. 
Various issues pertinent to system renewal are integrated into a common analytical 
framework to make for a more comprehensive and better formalised concept of 
obsolescence in military systems. 

This work is intended as a background study for later development of a com- 
prehensive policy to complement life-cycle management initiatives in the Singapore 
Ministry of Defence. It serves as a rudimentary framework that can be expanded 


later for actual use. 


E. IMPETUS 


The choice of this research topic was motivated by the following: 


e Singapore’s defence buildup has been underway for two decades and already 
system renewal decisions are due for some of the aging systems that were ac- 


quired initially. A framework that facilitates taking the long-term perspective 


in system renewal planning will help ensure that such optimal decisions are 
made in a systematic and timely manner. As the assets involved are consid- 


erable, improvements in system renewal planning is of significant importance. 


e Milestone 5 (upgrading/replacement) was introduced in U.S. DoD Instruc- 
tion 5000.2 in September 1987. No clear guidance has been provided on the 


approaches to be adopted for this new milestone. 


e Available literature generally either deals with replacements aimed at pre- 
empting failure (preventive maintenance) or reduces the problem to a purely 
economic one (commercial approach). There is a dearth of literature of rele- 


vance to military system renewal planning. 


F. SCOPE 


Replacement policies concerning routine preventive maintenance are not ex- 
amined in this thesis. The scope is limited to major systems or forces for which 
the considerable effort involved in identifying, tracking and analysing capability re- 
quirements, costs and measures of effectiveness (MOEs) are justifiable. At the heart 
of the framework to be presented is an analytical model, the workings of which 
will be illustrated under various assumptions based on a hypothetical weapon sys- 
tem. Attention is focused on single mission systems in changing environments. The 
treatment of threat evolution, force level planning and higher level optimisation 
are briefly discussed. Simple sensitivity analysis is illustrated but full and proper 


stochastic modelling and risk analysis lie outside the scope of this thesis. 


II. ISSUES AND LITERATURE SURVEY 


Systems are renewed as they become obsolete. Why? When is a system 
considered obsolete? What are the issues involved in planning for system renewal? 
These are the larger questions addressed in this chapter. We begin by explaining 
obsolescence and its consequences. Next, for a broad overview of the system renewal 
problem, various issues pertaining to the selection and scheduling of alternatives are 
discussed. The chapter ends with a survey of some of the research that have been 


carried out on system renewal. 


A. OBSOLESCENCE 


Obsolescence may be categorised as follows: 


e Economic: A system is considered economically obsolete when it becomes less 
cost-effective than available alternatives because of cost escalation attributed 
to physical deterioration, scarcity of spares, inflation, etc. An old truck that 
breaks down often, guzzles fuel that has become expensive, or has little spares 


support is economically obsolete. 


Technological: his form of obsolescence occurs when innovation makes avail- 
able a new, more cost-effective alternative to an incumbent system. Fighter 
planes with propellers, for example, became technologically obsolete when jet- 


propelled ones became available. 


Operational: Operational obsolescence arises wlien a system fails to meet its re- 


quirements because of deterioration in performance, enemy counter-measures 


or shifts in own strategic/tactical doctrine. Immobile defence systems, for 


instance, may become obsolete if it is known that they can be effectively en- 
gaged or bypassed by the enemy, or if strategic emphasis shifts from attrition 


to maneuver warfare. 


Obsolescence may be a compound effect of the above. A gun that is no longer 
accurate may also be uneconomical to maintain. 

The continued operation of an obsolete system is always a waste of money, 
but technological and operational obsolescence may have more serious consequences. 
Failure to exploit technological opportunities or remedy operational deficiencies may 
lead to military disasters. Technology is not exploited merely by keeping pace with 
research as was shown by the French experience in 1939. The French Army then 
had for some time developed prototypes of aircraft, armor and antitank weapons 
superior to those of the German Army. Their wait for still better ones to come along 
before introducing them into the inventory was a reason for their unpreparedness 
(Possony and Pournelle, 1970]. The Polish Army’s fatal cling to cavalry is a clear 
lesson in the consequence of operational obsolescence. 

Not all obsolescence is bad. The term “built-in obsolescence” is sometimes 
used for the practice of deliberately designing machines that deteriorate rapidly 
after a given life; the user being compensated for this by a lower purchase price. 


This can be seen, for example, in commercial batteries with different lives. 


B. STRATEGIC ISSUES 

Defence forces exist for a purpose, namely to deter aggression and, if deterence 
fails, to successfully wage war. The strategy to achieve this purpose should guide 
the formulation of system renewal plans. Potential adversaries should be identified 


and their intentions, war-proneness, capabilities, weaknesses and strategies assessed 


to arrive at a strategy that is congruent with foreign policy. Competitive strat- 
egy (Englund, 1987], for example, seeks a favorable cost exchange ratio by making 
obsolete major components of the enemy’s defence system. Stealth bombers and 
SDI are manifestations of this potentially destabilising strategy. Some references on 
strategic analysis and forecasting are provided in the bibliography. 

There are many facets of strategy that affects system renewal plans. Resources 
for modernisation may be channelled to boost readiness if conflict seems imminent. 
The danger of provoking or worsening an arms race must be weighed carefully if plans 
call for a significantly increased capability. The impact of a plan on the economy may 
also be significant. The B-2 program recommended by the USAF, for example, asks 
for a whopping seventy-five billion dollars. The impact of modest levels of defence 
spending on the economy is usually seen to be beneficial. Fredericksen’s [1989] 
analysis, for example, showed a favorable feedback relationship between defence 
spending and economic growth in Indonesia for the period 1964 to 1985. Economic 
issues and methodology in arms race analysis are surveyed in Leidy and Staiger 
[1985]. 

The link between strategy and procurement is historically lacking in the US 
(and Singapore). The chairman of a Defense Resource Management Study in 1979, 
for example, declared that there is broad agreement that the first P in PPBS is 
silent [Rice, 1979]. The Packard Commission in 1986 stated: 


... there is a great need for improvement in the way we think through and tie 
together our security objectives, what we spend to achieve them, and what we 
buy. The entire undertaking for our nation’s defence requires more and better 
long-range planning. This will involve concerted action by our professional 
military, the civilian leadership of the Department of Defense, the President, 
and the Congress. 
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Today, there is no rational system whereby the executive branch and the 
Congress reach coherent and enduring agreement on national military strategy, 
the forces to carry it out, and the funding that should be provided - in light of 
the overall economy and competing claims on national resources. The absence 
of such a system contributes substantially to the instability and uncertainty 
that plagues our defence program. These cause imbalances in our military 
forces and capabilities and increases the cost of procuring military equipment. 


This thesis seeks to improvise a link between strategy and force/resource plan- 
ning by incorporating strategic outlook into system renewal plans. The cost and 
effectiveness of alternatives of interest are projected into the future to make for 
more forward-looking plans. Figure 2.1 illustrates the links between strategic issues 
and system renewal planning. It is hoped that at some future time system renewal 
plans can be aggregated and linked to arms race models, strategic analysis models, 
econometric models, etc. to provide a higher level of optimisation. 

Technology pushes system renewal because it holds the key to one’s edge over 
the enemy. However, excessive emphasis on performance exacts an inordinate cost. 
Figures 2.2 and 2.3 [Gansler, 1989] show the cost histories and performance growth 
respectively of a series of U.S. fighter and attack aircraft. The performance assess- 
ment methodology used in Figure 2.3 is explained in Chapter 4. As can be seen, 
the difference in unit cost between high-end aircraft (e.g. F-15A, F-18) and low-end 
ones (e.g., F-16A, A-7A) had grown significantly through the years. Performance 
growth was however much more modest, indicating that cost-effectiveness of high- 
end aircraft had deteriorated relative to that of the low-end ones. Figure 2.3 suggests 
that the performance of the first operational units of a low-end aircraft lags that of 
high-end ones by only five years. It also implies that technological superiority over 
Soviet counterparts could have been maintained using lower-risk, lower-cost aircraft. 


The lower cost translates into more aircraft, which may matter more in terms of 
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military results. Expensive high-performance or multi-role assets also require in- 
tricate defences as they are too valuable to lose. Quality must also be traded off 
against quantity when selecting alternatives in a system renewal plan. This tradeoff 
can be accomplished within the model to be introduced in the following chapters by 
denoting quantity as a performance parameter. The impact of future innovation on 
the cost-effectiveness of alternatives of interest is projected so that optimal numbers 
of the right alternatives will be procured at the right time. 

The question of how far to push technology in a procurement program has 
major cost and effectiveness implications. Figure 2.4 illustrates how the technolog- 
ical frontier for a certain type of system may advance with R&D effort. Point A is 
typical of a target set to ensure the ability to counter probable threats over a long 
horizon. Due to the risks involved, it is likely in this case to end up somewhere inside 
the area bounded by the points A, B and C. The less ambitious and better plan 
would be to seek to quickly and cheaply arrive at say point D first, then improve 
incrementally to point A. Having such pre-planned product improvements (P*/) 
requires good forecasts of likely future upgrades. The model to be introduced facil- 
itates the selection and scheduling of such upgrades. The reader may be interested 
in the 1984 report by the U.S. Defense Science Board entitled “Improved Defense 
through Equipment Upgrades: The U.S. and Its Security Partners”. 

Requirement specifications largely determine the technological content of a 
program. These requirements may not be the result of careful study. In fact, they 
understandably tend to reflect the psychological comfort levels of those who have to 
deal with the exigencies of war. After all, it never hurts to be able to go the extra 
mile. Accordingly, the best of various candidate systems often becomes the “re- 
quirement”. The importance of specifying objective, mission-oriented requirements 


cannot be overstressed. 


Technological Economic Political Threat 
Environment Environment Environment Environment 


- = Present 
2 Operational 


Force 
















Arms Race 


Foreign Policy Threat/ Mission 


Force 
Formulation 


Structure 


Modelling Analysis 


Tactical 


Military Strategy Doctrine 


Formulation 

















Technological 
Forecasting 





Net Assessment 
( Deficiencies 
& Risks ) 


Ye OT & E, Insp., 











Economic 
Modelling & 
Forecasting 


Exercises & 
Wargames 







Force/ 
Renewal 
Planning 






Resource 
Planning 


Cost 
Models 











Prospective 
Future 
Alternatives 














Cost 
Models 


Functional Effectiveness Legal, Organizational 


Models 


Physical 
Models Models Budgetary, lemporal 


& Other Constaints 





Figure 2.1 : System Renewal and Related Issues 
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C. ECONOMIC ISSUES 


Cost-effectiveness of alternatives are dependent on future interest rate, foreign 
exchange rates, inflation rate, resale/salvage/scrap market condition and budget 
allocations. The U.S. Department of Defense publishes forecasts of inflation indices 
for various categories of military hardware which can be used in system renewal 
planning. DoD would of course like to estimate conservatively (high inflation) so 
as to get more money to play safe in case it actually is high. However, in times of 
high inflation, there is in fact considerable political pressure for lower estimates in 


order to suppress defence budget growth projections. The cumulative inflation from 
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1978 to 1981 was 45 percent, while the inflation estimated in the DoD budget was 
only 28 percent [Gansler, 1989]. Therefore one must exercise discretion in using the 
inflation indices. 

Assuming a common interest rate for all programs does not make for equitable 
comparison as alternatives with different life-spans and cost profiles are affected 
differently. One should try to forecast the rate as accurately as possible to avoid 
biased results. As economic factors such as interest rates, foreign exchange rates 
and market conditions are very difficult to estimate, sensitivity of solutions to their 
likely fluctuations should be tested in actual system renewal planning. 

Adoption of a longer term perspective in system renewal planning should be 
matched by the adoption of a multi-year budget. Procurement of needed quantities 
of a system is usually most economical at a certain rate. An optimal system renewal 
plan may involve humps in budgetary outlays. The minimum efficient rate of pro- 
duction for F-15 aircraft, for example, is 120 per year. However, DoD was allowed 
only 41 per year [Gansler, 1989]. This U.S. practice of stretching out programs to 
smooth out the humps tends to lead to gross inefficiencies in procurement. The 
floating of treasury bonds, for example, may be a better way to achieve a more 
efficient procurement program. For multi-year budgets to work, accurate cost esti- 
mates are necessary. Buying off-the-shelf items, adopting low-risk approaches (e.g., 
by emphasizing P*/) and designing to cost are collaborative means towards this 


end. 


D. MODELLING ISSUES 
Military system renewal planners have to grapple with problems of multi- 
facetedness, vagueness, uncertainty and subjectivity. Even if the effectiveness of a 


given alternative can be summarised by a single MOE, one has still to consider its 
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cost and probable service life. The system renewal problem has therefore at least 
three dimensions — effectiveness, cost and useful life. Often, service life 1s simply 
assumed and net cost is amortised into a rate over this period for evaluations. How- 
ever, advocates of an alternative will tend to assume a longer life than opponents. 
Instead of assuming when an alternative will become obsolete, the model to be pre- 
sented derives the service life as an output based on economic, technological and 
military (requirement) forecasts. This should make for more realistic and equitable 
evaluations. 

Risks due to uncertainty can be assessed by testing the sensitivity of solutions 
to possibilities learned from experience. Dominance or robustness of a solution 
would provide some assurance of its validity under a range of expected scenarios. 
The quality of experience is of critical importance for such risk assessment. As such, 
qualified experts should be consulted in the formulation of scenarios. 

Besides multi-facetedness and uncertainty, a planner has to contend with prob- 
lems of vagueness. Strategic analysis may highlight the need to heighten readiness in 
the coming years but such qualitative goals do not readily lend themselves to math- 
ematical modelling. Even if readiness can be precisely quantified by a parameter, 
there may not be consensus concerning what constitutes the proper level of readi- 
ness. Fuzzy set and logic theory offer means to bridge the gap between quantitative 
and qualitative reasoning, thus linking strategy formally to procurement planning. 
The efforts of J. Dockery [1985] is an example of work done in this area. Karprzyk 
and Orlovski [1987] and Zimmerman, et al., [1984] also discuss new paradigms more 
suited to situations where vagueness and uncertainty are clearly present. The prob- 
lem of vagueness 1s not considered in this thesis but may be addressed in future 


embellishments. 


Although the exercise of human judgment is essential in system renewal plan- 
ning, the subjectivity involved must be managed well to reduce the possibility of bias 
(conscious or otherwise). The problems encountered in assessing utility functions, 
for example, were surveyed by Hershey et al., [1982] and Farquhar [1984]. It should, 
however, be noted that intuitive (subjective) judgment of experts may in certain 
cases prove superior to analytical studies [Hammond, 1987]. Human behavioral as- 
pects are important in a real application, but they are not considered in this thesis 
because the cases studied are hypothetical ones meant for illustrating the model 
that was developed. Utility functions, for example, were assumed to be explicitly 
defined even though such functions are difficult to assess accurately in practice. 

The horizon of optimisation is an important issue in modelling system renewal. 
Very often, this horizon is either fixed at some arbitrary value or taken to be in- 
finity. Both practices suffer from serious flaws. Choosing horizons that are a small 
non-integral number of life-spans would produce incorrect results. This occurs, for 
instance, if the horizon for a system with a typical lifespan of ten years is set to, 
say, twenty-five years. Assuming an infinite horizon is clearly unrealistic. There are 
procedures that allow one to find the horizon beyond which the decision for the first 
change remains optimal. However, these so-called “planning horizon procedures” 
are essentially also based on the infinite horizon assumption. This thesis takes a 
different approach, namely the finite-stage approach (from dynamic programming), 
in which the optimisation horizon is a variable output given by the years spanned 
by a sequence of systems. 

Modelling the effectiveness of a military system is sometimes confused with 
predicting the outcome of its use in battle. The two are related but not synonimous. 
As an example, let us suppose that a certain armour unit has the MOE structure 


shown in Figure 2.5. The outcome in a given battle scenario may be described 
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by the number of own and enemy losses. Neither of these numbers expresses fully 
the effectiveness of the force. An abstract figure of merit would have to be com- 
puted from these numbers based on judgment concerning what constitutes success. 
Therefore, one might say that battle outcomes are predicted by combat modelling 
whereas effectiveness is evaluated by “combat value” modelling. The abstract figure 
of merit may be said to be the MOE, and performance parameters that determine 
this value may be said to be the measures of performance (MOPs). This distinction 
is, however, an artificial one as the MOPs may collectively be considered to be a 
vector-valued MOE. Accordingly, MOE and MOP in this thesis are both taken to 
mean the same thing — key performance parameters that determine the effective- 
ness of a system. In the example shown in Figure 2.5, the MOEs would be range, 
accuracy and lethality if these parameters are the ones that have major effects on 


effectiveness. 


E. ORGANISATIONAL ISSUES 

System renewal plans may be developed for many different levels, ranging from 
a single equipment to an entire task force. Unless they are coordinated such that the 
objectives of upper echelons are considered by lower echelon planners, sub-optimal 
or even counterproductive plans are likely to be made. Most resource allocation 
models do not consider the organisation as a factor in the decision process; that 
is, they assume a monolithic decision-making entity. The same decision-makers are 
assumed whether the subject of planning is a single artillery gun or an infantry divi- 
sion. Ruefli [1974] surveyed several economic and management science models that 
provide remedy for this unrealistic assumption by explicitly reflecting organisational 
structure in the model. The Dantzig-Wolfe model [Dantzig and Wolfe, 1960] is the 


best known of the classical models. In more recent years, models such as the one by 
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Figure 2.5 : Hierarchy of MOEs 


Source : Modified version of the diagram from 
“Notes on MOEs” by E.B.Rockower, Aug. 1985. 


Cassidy et al., [1971] that are based on behavioral theories have emerged. Organ- 
isational issues are not considered in this thesis as the model is still a conceptual 


one that is not ready for organisation-wide implementation. 


F. LITERATURE SURVEY 
No literature suitable for military system renewal planners could be found as 
all those surveyed were commercially oriented ones that subjugated performance 


considerations under cost. Most of them consider only replacement alternatives. 
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The earliest of these were based on simple calculations and differential and integral 
calculus. Examples are Preinrich’s Constant Chain Model [1940] and Terborgh’s 
MAPI system [1949]. In more recent years, methods such as dynamic program- 
ming, Markov chains and Pontryagin’s maximum principle have become more pop- 
ular. The first such model was introduced by R. Bellman [1955] using dynamic 
programming. Howard [1960] used a policy-iteration method to solve multi-stage 
decision problems in connection with Markov processes. Numerous articles conern- 
ing stochastic models were published in the 1960’s. A survey of the treatment of 
stochastically failing equipment was compiled by Jorgenson et al., [1967]. Naslund 
[1966] was the first to apply the maximum principle to the replacement problem. 
Thomson [1968], Rapp [1969], Sethi [1970] and many others used these so-called 
“optimal control” theory models that uses the maximum principle. There are nu- 
merous variants in replacement models, the assumptions of which were grouped by 
Luxhoj and Jones [1986] in a general framework. Some of the works not referenced 
in this survey may be found in the bibliography. 

Some of the more noteworthy recent publications that explicitly considered 


technological advances are as follows: 


1. R.C. Stapleton et al., [1972] examined the effects of different forms of tech- 
nological change (change in purchase price, resale value, operating costs and 
rate of deterioration) on costs using dynamic programming. This work pro- 
vided the author of this thesis with basic ideas on how to model the impact 


of technological progress on cost. 


2. S. Chand and S. Sethi [1979 and 1982] applied the planning horizon procedure 
developed by H.M. Wagner and T.M. Whitis [1958] for the dynamic lot size 


model to the replacement problem. It allows one to obtain the first replacement 
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time for an infinite planning horizon based on the optimal policy for a finite 
horizon. Technological advances is forecasted for this equivalent finite horizon. 
The 1982 paper, which extends the first to the case of multiple alternatives, 
treated upgrades within the same framework as replacement alternatives. This 
idea led to the adoption of a common framework for all alternatives in this 


thesis. 


. T. Goldstein et al., [1988] examined the influence of,a forthcoming new and 
better machine (based on an expected technological breakthrough) on the 
system renewal decision. Each year, if the new machine does not appear, a 
decision is made to either continue with the incumbent and wait for the bet- 
ter machine, or buy a new version of the incumbent. The model determines 
the optimal age for replacing the existing machine. This model ignores evolu- 
tionary changes that characterises technological innovation and thus was not 


adopted in this thesis. 


. Y. Kusaka [1985] showed a method for deciding whether to keep or replace 
the existing machine at the present time without determining the sequence of 
subsequent replacement times. A joint paper with H. Suzuki [1988] expanded 
on the earlier evaluation system by introducing the concept of “control limit” 
policy. Kusaka’s method was not adopted in this thesis as the provision of 


estimates concerning future renewals is deemed to be a useful feature. 


nS 


Ill. THE MODEL 


A. INTRODUCTION 

As evident in the previous chapter, many complex issues impinge upon de- 
cisions concerning systems renewal planning. Literature searches have failed to 
uncover any model suitable for system renewal planning in the military. The dy- 
namic system renewal planning (DSRPM) model to be presented in this chapter 
is a first-cut attempt to weave a unifying framework spanning some of the more 
quantitative issues involved in systems renewal planning. The overall approach is 


outlined first, followed by theoretical underpinnings of the model. 


B. APPROACH 

The adopted approach is to use a versatile methodology, namely dynamic pro- 
gramming augmented by the Lagrangian relaxation technique, to explicitly model 
the effects of aging, technological advancement, and changes in capability require- 
ments. The methodology operates on an acyclic directed network of the kind shown 
in Figure 3.1. This network representation of the system renewal process is unique 
to this thesis. As this network lies at the heart of the model, a brief description is 


in order before we move on to the other parts of the model. 
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Note 1: Nodes denote possible transitions while arcs 
represent possible alternatives and system life-spans. 
The dotted arc from node 3 to 7, for instance, stands 
forthe possibility of acquiring alternative 2 in1990 and 
using it as System B until 1992. 

















Note 2: This diagram is just an illustration. In real 
applications, the number of arcs, nodes and years 
spanned would be higher. The years spanned by 
the network also need not be contiguous. 


Note 3: System A and Alternative 0 are 
synonymous with the present system. 
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Figure 3.1 : Network Representation of System Renewal Process 
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The network represents the system renewal process involving the selection of a 
succession of systems to fulfill some capability requirements. The root node 1 repre- 
sents the present system (System A) in the current year. Each of the arcs emanating 
from it represent a possible residual life-span of this system. These arcs terminate 
on a column of nodes representing possible points in time for transition to the next 
system (System B). Transitions may be replacements or upgrades. Similarly, arcs 
extending from these nodes to an adjacent column of nodes represent the possible 
years for transition from System B to System C. The number of arcs connecting each 
pair of nodes corresponds to the number of alternatives that may be selected for 
System B; likewise for System C. Thus if cumulative cost and effectiveness of using 
a given alternative over a given life-span are assigned as attributes to the associated 
arc, optimisation for the three systems reduces to one of finding the “best” path 
through the network. 

Finding the “best” path constitutes solving a three stage dynamic program- 
ming problem. If cost is set to a very large number whenever an alternative ages 
such that one or more performance parameters at the end of a given life span fails 
to meet stipulated requirements, then dynamic programming yields the solution 
that meets the requirements at minimum cost. However, if one seeks the optimal 
trade-off between cost and effectiveness, then the problem becomes a multicriteria 
one. This problem can be solved using the Lagrangian relaxation technique by for- 
mulating this problem as one of maximising effectiveness given a certain budgetary 
limit on cost. Details of this technique are given in section D of this chapter. The 
“best” path can also be solved using linear programming formulation of the network 
flow problem. However, the dynamic programming formulation has more intuitive 
appeal and lends itself easily to incorporation of lead-time, discounting and other 


considerations. 
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The aging process is modeled by discounting performance and cost character- 
istics over time, based on empirical parametric relationships. What these charac- 
teristics will be for a given alternative acquired new in a certain year is projected 
based on technological trends. Forecasts of inflation and foreign exchange rates are 
also factored into the projected acquisition, maintenance and other costs. All out- 
lays are converted to current year dollars based on forecasted interest rates. The 
projected rates may vary from year to year but a steady trend is assumed beyond an 
arbitrary forecast horizon. Performance requirement levels may also be stipulated 
for each year based on analysis of politico-military trends. A unique set of life-cycle 
cost and effectiveness measures is thus associated with the network for each scenario 
expressing a given set of expectations concerning relevant future trends. The mini- 
mum cost solution and efficiency boundary for each scenario can then be obtained 
by applying dynamic programming and Lagrangian relaxation to the corresponding 
network. The efficiency boundary expresses the upper bound on effectiveness over 
the entire spectrum of life-cycle cost spanned in a given scenario (Figure 3.2). The 
point on this boundary that yields the highest effectiveness-to-cost ratio is the so- 
lution that gives the most bang for the buck. As the boundary is piecewise linear, 
all points along a line segment may qualify if extrapolation brings the line through 
the origin. Otherwise, a Lagrangian solution point at the intersection of two line 
segments will qualify. 

The reader may wonder why a three stage process is adopted in this model. 
The model itself does not restrict the user to any number. Rather, the matter 1s 
one of judgment concerning the credibility of long-term forecasts. Although it 1s 
desirable to optimise over as long a stream of systems as possible, it is also futile 
to try to visualise the distant future. A two-stage model comprising the present 


system and its immediate successor is satisfactory as far as a reasonable optimisation 
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horizon goes. A three-stage model would be preferable as it allows realistic bounds 
to be developed for the immediate successor’s economic-technological life. The cost 
of competing alternatives can therefore be more equitably compared since their 
amortisation periods may be more accurately estimated. Infinite horizon models 
meant for static environment problems are clearly not suitable for our purposes. 
Accordingly, assumptions of monotonicity, contraction, boundedness, etc., which 
are required for asymptotic stability of such models [Henig 1978, Beckmann 1968, 
Boudarel 1971] need not be made. 

It may be noted that a long lead time often precedes transitions, be they 
upgradings or replacements. Lead time of an alternative is defined to mean the 
minimum number of years by which the initiation of acquisition must precede the 
achievement of initial operating capability (IOC). By setting to infinity the cost of 
alternatives that cannot be introduced because of insufficient preceding lead time, 
we can avoid generating solutions contradicting lead time requirements. This check 
is made not only for transitions from the present system to its immediate successor, 
but also for subsequent ones. The latter requires the assumption that initiation 
of a transition never begins before the preceding one even takes place. A path in 
the network is thus deemed infeasible if there is any arc representing a system life 
shorter than the lead time required by its successor. 

While no one can reliably predict the future, one ignores future possibilities 
and trends only to one’s disadvantage. Taking calculated risks after examining 
a broad range of scenarios is the only realistic approach to handling uncertainties 
about the future. For this purpose, the flexibility of the model presented is a valuable 
asset. Long-range projections can be varied more than near-term ones to account 


for greater uncertainty in these projections. 
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It should be emphasized that the focus of analysis is on the timing and choice 
of alternatives for the first transition, the optimality of which is affected less by 
inaccuracies in long-range forecasts than short-term ones. By updating the scenarios 
annually for changes in the perception of qualified experts, rolling plans can be 
developed for system renewal decisions to be made optimally and in a timely fashion. 


For this reason, the model is named the Dynamic System Renewal Planning Model 


(DSRPM). 
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Figure 3.2 : Domination of Solution Space by Efficiency Boundary 
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ASSUMPTIONS 


For the purposes of this thesis, the following assumptions have been used: 


A steady state exists beyond an arbitrary forecast horizon where all projections 
are either constant, linear or harmonic. This assumption expresses the lack of 


detailed knowledge about the remote future. 


Cost and performance parameters can be measured deterministically and util- 
ity functions can be obtained explicitly. (Note: This assumption would be 
relaxed in future embellishments incorporating stochastic models and implicit 


approaches to utility modeling.) 


An upgraded system is a new system with acquisition cost equal to the direct 
upgrading costs plus the salvage or resale value of the original system in the 
previous year. A new acquisition with no precedent is simply one in which 
the incumbent system is a null or dummy system. With this assumption, 
upgrading, replacement and new acquisition alternatives are integrated within 


a common framework. 


Acquisitions (upgrading or replacement) and their associated outlays occur 
at the beginning of a transition year. The old system is disposed of at the 
end of the year. This definitive assumption provides for proper costing of 
the metamorphic process by which a new system achieves initial operating 
capability (IOC) while the old is prepared for disposal. Disposal does not 
really occur in an upgrading, but the new capability is still assumed to be 


achieved at the end of the transition year. 


No transition is initiated before the preceding one takes place, as was explained 


in the previous section. 
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One should note the small number and unrestrictive nature of the assumptions 
made in the model. This is a direct consequence of the versatility of the methodology 
employed. There are other assumptions needed to model particular alternatives and 
scenarios such as those pertaining to utility functions, nature and impact of trends, 


etc., but these are not generic to the model discussed. 


D. METHODOLOGY 
Dynamic programming and Lagrangian relaxation are the cornerstones of the 
model’s integrity and the source of much of its power. This section explains these 
key techniques as applied to finite discrete decision processes. As methodologies 
concerning forecasting, cost estimation and other components of the model may 
vary according to the peculiarities of the application, these will not be discussed. 
The reader is referred to the bibliography for materials on these aspects of the model. 
1. Dynamic Programming 
Borrowing the notation used by the pioneering developer of dynamic pro- 
gramming, Richard E. Bellman, [Bellman and Kalaba, 1965], we begin by examining 
the important notion of policy in the context of system renewal. If a decision d; in 
the 2" stage results in a transition from a state y; to a state y;41, then an n-stage 


process may be represented as 


Ole T' (Yo, do) 
a T(y1, 41) 
Yn = V(Ga-idee) (3.1) 


In the 3-stage system renewal process depicted in Figure 3.1, for instance, 


the process of operating the present system for a year, then acquiring alternative 
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1 to use for two years before acquiring anew another alternative 1 for a further 
two years use would be represented by a path P\,37,2 of unbroken line arcs passing 
through nodes 1, 3, 7 and 12. 

Suppose we are concerned with maximizing a prescribed scalar function 


of the states and the decisions: 


Ryo. ¥15 Y20° °° Uns dy, d,,°-+,d,_}) (3.2) 


This so-called criterion function, in the case of a system renewal process, may be 
the cumulative expected return accrued from acquiring and operating a particular 
succession of alternatives, i.e., traversing a given path in the network. Without an 
optimal rule as to how decisions should be made, the only way to obtain the max- 
imum of /? and its associated path is by exhaustive enumeration along all possible 


paths. Effort is much reduced with the use of rules in the form 


dy, = dy(¥o. 15° °°) Ye: do, dy,:--,dx_1) (3.3) 


These functions are called policy functions, or simply policies. A policy that maxi- 
mizes the criterion function f is called an optimal policy. 
Let us consider, in particular, policies which depend only on the current 


state, k, that is policies of the form 


dy = dk (yx) (3.4) 


This additional simplification is possible only when the structure of the criterion 
function R has the important property of divorcing the past from the future. I5x- 


amples of criterion functions with this “separable criterion” property include 
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2 Cyd) (3.5) 


and I r(y;, 4;) (3.6) 
The first example is of die most common type of criterion. The usual criterion of 
cumulative cost or cumulative returns used in system renewal is of this variety. The 
reliability of series systems would perhaps come to mind in the case of the second 
one. An example of a criterion function that does not have the separability property 


is that of an overall ratio, that is, 


E(yi, d:) 


Me 
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(3.7) 
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C'(yi, di) 
The use of overall cost-effectiveness ratio as a criterion in the renewal process would 
thus require decisions to be based not only on the current state but also on the path 
leading to that state. 

When the criterion function used has this vital separability property, the 


optimal policy is characterized very simply by the following principle of optimality: 


(Qenvise sums tmaie te 1S optimal > (Yo: Y15°'*y Ym) 1S optimal (3.8) 


This intuitive result can be proved by contradiction. Suppose that for the (m + 1) 
stage process, (Yo, 41,°°';Ym) 1s not optimal. Specifically say (29,2),°°:,%m) is 
optimal. Then (29, 21, --+,; 2mstimeiy 7 oes better tham (yoriy°*:., vente 
---, Yn). But this contradicts the original preinise that (yo, y1,°°*, Ym) Ym419'* 5 Yn) 
is optimal. Therefore, expression 3.8 must be true. Using the principle, one obtains 


for the examples given by expressions 3.5 and 3.6 the following optimal policies: 
Skai(T(yn,dk)) = max [c(ye, de) + Se (ye)] (3.9) 
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Seas(T (ye, dk)) = max [c(ye, dk) fi (ye) (3.10) 


where f;(y;,d;) is the criterion function in state 7 obtained using an optimal policy. 

Dynamic programming hinges around the formulation of optimal policies 
to enable solutions to be found by systematic recursion (Beckmann, 1968]. A solution 
may be given in terms of the sequence of criterion functions { f;(y;)} or the sequence 
of policy functions {d,(y,)} . Each determines the other, although one should note 
that there may be several optimal policies yielding the same sequence of criterion 
functions. A policy may also be satisfied by one or more solutions. Only when the 
solution is unique is the policy a necessary and sufficient condition for optimality. 
Otherwise, it is necessary to compare the solutions using additional considerations 
such as risk, expected life-span, etc., to define the optimal solution. 

Being a flexible and potent tool, dynamic programming has many vari- 
ations. In stochastic dynamic programming, the next state is chosen based on 
conditional transition probabilities and the criterion function is given as a random 
variable, usually the mean of possible rewards [Ross, 1983]. Variations with vector- 
valued criterion functions have been developed [Henig, 1978], but these multicriteria 
dynamic programming models are still not as widely applied as the stochastic vari- 
ants. Although these multicriteria models could be used in system renewal problems 
to obtain dominating (Pareto-efficient) solution sets, it is simpler, more useful, nat- 
ural and elegant to use the Lagrangian relaxation technique in such problems, as 
shall be shown in the following section. 

2. Lagrangian Relaxation Technique 

Two conceptual approaches used in systems analysis are either to fix 

effectiveness and minimize cost or fix the budget and maximize effectiveness [Fisher, 


1970]. Suppose we choose the latter and formulate the system renewal problem as 
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ymax 1 Del (y;,d (at) 


subject to 


ror yi, d;) < Budget (al) 


where E(y,,d;) and C(y;,d;) are the life-cycle effectiveness and cost respectively 
of making a decision d; in state y;. Units for the former are utiles, while that of 
costs are dollars. In the system renewal process shown in Figure 3.1, for example, 
if the current state is represented by node 3, and the life-cycle effectiveness and 
cost of acquiring alternative 1 and operating it for two years are 2,100 utiles and 
$7m respectively, then /(3,2) = 2,100 utiles and C(3,2) = $7m. The constraint 
in expression 3.12 may be relaxed by moving it to the objective function with the 


inclusion of a Lagrangian multiplier 4, thus 


jmax »» E(y;, d;) — 2 )- a1" C(y,,d;) — Budget | (3.13) 


Noting that expressions 3.5 and 3.13 are of the same form, we know that the optimal 
policy for maximizing this criterion function is given by the principle of optimality 
as: 
Sasi (T (yas de) = max {[E (ye, di) — AC(yes de)} + Se(ye)} (3.14) 
Therefore, for any fixed A, conventional dynamic programming may be used to 
obtain the optimal solution to expression 3.13. Lagrangian relaxation is tlius a 
technique for converting constrained optimisation problems into unconstrained max- 
imization (or minimization) problems [Everett, 1963]. It may also be interpreted as 
a technique to solve the dual of the original problem given by Equations 3.11 and 
3.12, as will be apparent in the following discussion. 
To consider the effect of varying A, let us define the Lagrangian function 


L(A) to be equal to the objective function given in expression 3.13. That is, 


32 


HEN) = max >> E(y;,4;) — A So Cluid (y;,d ) ~ Budget] | (3771S) 


Oxy dn 0 
This piecewise linear function is illustrated in Figure 3.3. 
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Figure 3.3: Lagrangian Function L(i) 
When A = 0, the budgetary constraint is ignored and we obtain the 
maximum effectiveness solution to the (n + 1) system renewal problem given in 


expression 3.13. Denoting the cumulative life-cycle effectiveness and cost of this 


solution as LC Eg and LCC» respectively, we note that L(0) = LCL and that the 
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slope at this vertical axis intercept is —(LCCo — Budget). For marginal increases in 


A beyond zero, 


L(A) = LCE — (LCCy — Budget) (3.16) 


which is the equation of the leftmost line segment. When X reaches a threshold 
value, the solution changes to one with smaller cumulative life-cycle effectiveness 
and cost. The next line segment will thus have a lower intercept value and an 
increased slope. This therefore implies that the Lagrangian function will always be 
convex for system renewal problems formulated as in expression 3.13. 

The Lagrangian multiplier \ has a physical interpretation, namely that 
of the relative value of cost and effectiveness, i.e., shadow or dual price. Its unit is 
utiles per dollar. As 4 is increased from zero to infinity, the Lagrangian function 
steps through the entire solution space, from the maximum effectiveness solution to 
the minimum cost solution. The latter corresponds to the rightmost line segment 
which extends to infinity. What then is the optimal value of 1? To see this, we first 
note that for \ > 0, L(A) is an upper bound on the optimal value F* of the original 
objective function given in expression (3.11). This is because minL() is the dual 


of the original problem. The minima of L() therefore yields F*. That is, 
ae DOM 10 Nai (sole 


This is akin to seeking the constrained stationary points in analog applications of 
Lagrange’s method to differentiable functions. 

The minima of L(A) do not necessarily exist. In a commonly observed 
case, LCCo 1s higher than the budget and the slope at the vertical intercept 1s 
negative. As D(X) is convex, A* is always finite and non-zero in such cases. If, 


however, LCCpo is lower than the budget, that is to say the maximuin effectiveness 
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solution is affordable, then the slope at the intercept will be positive and we obtain 
A* = 0. Should LCC be higher than the budget and there are no affordable 
solutions, then there will be no line segments of non-negative slopes, no minima and 
hence no feasible solution. 

The value of \* can be found most expeditiously and accurately by the 
supporting plane method illustrated in Figure 3.4. Given the upper and lower 
bounds on A*, the method involves the extrapolation of the line segments corre- 
sponding to these bounds to obtain the intersection value \;. The intercept and 
slope of these line segments are given by the cumulative life-cycle effectiveness and 
cost respectively of the solution obtained by solving expression (3.13) for the bound- 
ing values of 4. The lower or upper bound is replaced by 4; if the slope of the line 
segment corresponding to 4; if the slope of the line segment corresponding to 4; is 
negative or non-negative respectively. This process is iterated until \; equals one of 
the bounds. This equality occurs because L(,) is piecewise linear and convex, which 
enables \* to be obtained precisely. The starting lower bound 1s set to be zero. Any 
suitably large value of 4 that yields a line segment with non-negative slope may be 


used as the initial upper bound 4,,,,, where 
a NOMONETLOCGS eG BOE a. Onin (3.18) 


Experience of the author with similar problems shows this supporting plane method 
to be faster and more accurate than the commonly used bisection method, and A* 
may be found in only a few iterations even for very large networks with thousands 
of nodes and arcs. 

For each scenario, using the methodology described, we are able to obtain 
the solution points marking the efficiency boundary shown in Figure 3.2. This 


is begun by artificially setting a sufficiently high budget to obtain the maximum 
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L(A) 
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Figure 3.4: Iterations for Derivation of \* 


effectiveness solution, which is marked 1 in Figure 3.5. The corresponding value 
of * is obviously zero. ‘The budget is then set to just under the value of the 
cumulative life-cycle cost C, of the maximum effectiveness solution. This produces 
the next solution point, marked 2 in the diagram, and the process is repeated until 
we reach the minimum cost solution, which is the sixth solution. No solutions will be 
obtained for budgets less than Cg. The optimal Lagrangian multiplier 4* obtained in 
each case determines the slope of the line segment to the right of the corresponding 
solution point. It may be noted that the efficiency boundary consists of bounding 


hyperplanes which envelope the solution space but do not pass through all the 
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Pareto-efficient points. There are, for instance, three such points denoted by circles 
for budgets between C3 and C, which will not be discovered by the Lagrangian 
relaxation method. The basic cause of these inaccessible regions (called gaps) is 
non-concavity in the relationship between the optimal effectiveness and cost. These 
gaps will not exist if the law of diminishing marginal utility holds strictly. H. 
Nakayama, et al., [Nakayama, 1975], has a good survey of methods to deal with this 
characteristic of the Lagrangian relaxation method. The existence of these gaps does 
not create a problem as far as the application of the method in systems analysis is 
concerned. This is because the solution offering a maximum overall effectivenss-to- 
cost ratio is always a Lagrangian solution point. Also, a straightforward application 
of dynamic programming will produce the minimum cost solution even if it does not 
lie on the efficiency boundary. Note that as system lives are modeled explicitly for 
economic, technological and politico-military changes, the net cost of alternatives 
can be amortized over the proper period to reveal the true minimum cost solution 


as defined by amortized annual cost, as will be amplified in the following paragraph. 
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Figure 3.5 : Derivation of Efficiency Boundary by 
Lagrangian Relaxation Technique 


We have seen how characterization of the systems renewal problem as a 
constrained optimization problem (in which one seeks to maximize effectiveness for a 
given budget) enables us to obtain the efficiency boundary corresponding to a given 
scenario. A band of these boundaries is found for each set of scenarios investigated. 
If probabilities are attached to scenarios, our “expected” efficiency boundary and 
its associated optimal solutions can be obtained. Furthermore, bearing in mind that 
the cost of an arc in the network is set to infinity if there is insufficient lead time 
preceding it or if any performance parameter deteriorates below its requirement level, 
we know that any solution with finite cost is feasible from the effectiveness point of 
view. The minimum cost, in current year dollars, for the various feasible solution 
paths, as well as the associated period of use, can then be translated into amortized 
costs flows for given interest rate assumptions. On this basis, we can find the solution 
that allows us to minimize equivalent annual cost while meeting a given set of lead 
time and performance requirements. Therefore, both the fixed budget and the fixed 
effectiveness approaches are embodied within the model, delivering the maximum 
overall effectiveness-to-cost ratio solution and the minimum cost satisficing solution 
respectively. Furthermore, it should be noted that A* has the interpretation of 
shadow or dual prices in economic analysis. ‘Therefore, efficiency boundaries of 
systems with similar MOEs can be meaningfully compared. Also, choice of system 
mix for various capabilities may be guided by the fact that shadow prices should 
be roughly the same for all the systems in an optimal mix. The combined use 
of dynamic programming and the Lagrangian relaxation technique thus provide us 
with a conceptually simple, yet potent and useful methodology that is completely 


congruent with and derivable from the tenets of systerns/economic analysis. 
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IV. IMPLEMENTATION 


A. INTRODUCTION 

The dynamic systems renewal planning model (DSRPM) presented in the pre- 
vious chapter provides a formal basis by which diverse economic, technological and 
politico-military concerns may be integrated within a common framework (Figure 
4.1) to make for more rational management of defense assets. Credible and valid 
results, however, can be obtained only with the active participation of relevant logis- 
tics, engineering, finance, intelligence and operational agencies in the development 
of proper effectiveness and cost sub-models as well as realistic scenarios. High-level 
management support is thus essential for the successful implementation of the model. 
To foster a better understanding of the model and to demonstrate its plausibility 
and usefulness, a small prototype was developed based on a simple hypothetical 
weapon system. This prototype is documented in the next chapter along with case- 
studies. Although the addition of selected successful applications would give greater 
credence to the model, the limited time available for work on this thesis precluded 
full-scale implementation for any real systems. Nevertheless, effort shall be made 
to outline conceptually how implementation is envisaged. Only aspects of cost. es- 
timation, effectiveness modelling and requirement formulation peculiar to DSRPM 
is discussed. For a full discussion, the reader is referred to relevant materials in the 
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Figure 4.1: System Renewal Planning Framework 
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B. CONCEPTUAL FRAMEWORK 

Implementation of the model is best explained in terms of its constituent el- 
ements and their relationships to the processes involved in defense planning. As 
depicted in Figure 4.1, DSRPM has five elements, all of which are part of the 
force/resource planning process. A model with these five elements is developed 
for each entity of force planning, which may be weapon systems, task forces, bat- 
tle groups, etc. In renewal planning parlance, all these entities are simply called 
systems, and candidate systems are called alternatives. 

The cost estimation element projects the life-cycle costs for both the present 
system as well as future alternatives. The effectiveness modelling element concerns 
the appraisal of the value of a system as a function of its measures of effectiveness. 
The deterioration modelling element accounts for the change in effectiveness and 
cost over the life of a system caused by wear-out, operational attrition and enemy 
response. The other two elements are self-explanatory. All of these elements and 
their implementation are elaborated in sections C, D, and E. Only deterioration 
modelling and the selection and scheduling of programs are organic to the DSRPM. 
The other elements are on-going force/resource planning processes adapted with 
future projections for use in the DSRPM. 

Formulation of future alternatives is a part of the force/resource planning pro- 
cess that lies wholly outside the province of the DSRPM. It involves definition of 
the physical and functional forms of alternatives identified in the strategy formu- 
lation process and is shaped by technological outlook. Alternatives encompass the 


following: 
1. Conversion of incumbent systems to serve different roles 


2. Upgrading of incumbent systems to serve expanded or higher roles 
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3. Modernisation of incumbent systems to better serve their assigned roles 


4. Replacement of incumbent systems with newer and probably better alterna- 


tives 


5. Activation of reserve systems for active deployment (incumbent systems are 


the ones in dormant or preserved status) 
6. Reduction of number of incumbent systems through disposal/sale 
7. Retirement of incumbent systems into the reserves 


8. Acquisition of new (no precedent) systems to meet qualitative and/or quanti- 


tative deficiencies or new (threat) requirements. 


No present system exists, i.e., system A in the network shown in Figure 3.1 


does not exist, in the case of new acquisitions. 


C. EFFECTIVENESS MODELING 

The effectiveness modeling adopted in the DSRPM is completely congruent 
with the way military systems are selected and managed. Each system has one or 
more missions that represent the raison d’etre of that system. These missions are 
usually phrased in terms that are too broad for concrete guidance. What it takes to 
carry out the mission(s) are stipulated clearly as required operational (capabilities) 
requirements. These are usually lower bounds set on key performance parameters 
such as speed, payload, survivability, operational range, etc. The effectiveness of 
candidate systems is evaluated based on these key parameters and some criterion. 
The decision to choose a particular system need not be based exclusively on objec- 


tive parameters, and the criterion used is not necessarily a formal one. The views 
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of qualified experts and decision makers are also important. In fact, the politi- 
cal process leading to the decision may even be colored by elements that detract 
from the military contribution to national security, such as amicable relations with 
another state. However, assuming rational arguments hold sway, an effectiveness 
model representing the embodiment of the values, preferences and possibly biases 
of the decision makers, tempered with experience gained from performance audits, 
can be constructed. Justification papers, test and evaluation plans, trial reports, 
etc. are good sources of information needed to construct an effectiveness model. 
Combat models and wargames also provide valuable input. The aim of effectiveness 
modelling for this model is to describe value judgements of the decision-makers and 
not to predict battle outcomes (although predictions do tend to affect judgement). 
Therefore, the same effectiveness model used for evaluation of alternatives in the 
acquisition phase, and possibly altered by performance audits during deployment, 
is incorporated in the DSRPM. The only difference in the DSRPM< case is that 
prospective alternatives are evaluated alongside the incumbent, and performance 
parameters are projected both for future generations of each alternative as well as 
through the years of use after an alternative is acquired. Figure 4.2 shows how the 
j*® performance parameter of an alternative as acquired in year y may change over 
the years of usage tf. 

As an example, let us suppose that a weapon system was evaluated based on 


an effectiveness model with the following additive structure: 
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where 
e TE, is the total effectiveness of alternative a. 


e W; is the weight of the j*" performance parameter expressing its importance 


relative to the other parameters. 
e X;, is the value of the j* parameter for alternative a. 
e U; {Xj,.} is the utility function that transforms X;, into a utility value. 


e X, is the minimum level required of the j** parameter (operational require- 


ment ). 


In this case, the effectiveness sub-model used in the application of DSRPM to the 


system will be: 


T Bary(t) = J W5U; {Xjarydars(t)} (4.2) 


j=l 
where 


e TE.) is the total effectiveness of alternative a’ (not necessarily the same as 


a) if acquired new in year y during the t*” year of deployment. 


e X;,’, is the value of the j performance parameter for alternative a’ if ac- 


quired in year y. 


e d, ;(t) is the deterioration function for the j* parameter of alternative a’ in 


the t* year of deployment. 


The alternatives a and a’ need not be the same. The latter typically comprises 
not only the replacement alternative, but also upgrading options. Note that the 


structure of the effectiveness model, weights and utility functions are retained in 


the DSRPM. TE,.,(t) is the effectiveness associated with all arcs representing al- 
ternative a’, leading out of nodes of year y, with life-spans of ¢ years in the network 
shown in Figure 3.1. If Xarj;(t) < Xjy41 for any j, i.e., if any parameter fails to 
meet the requirement level stipulated for the year of operation (y+t), then the cost 
of the corresponding arc is set to infinity. 

The effectiveness model used here is similar to that of the TASCOM method- 
ology developed by The Analytic Sciences Corporation for OSD over several years 
in the late 1970’s and the early 1980’s. The TASCOM model provides a first-order 
indication of an aircraft performance potential by computing a score based on at- 
tributes such as range, speed, maneuverability, payload, basing mode, target acqui- 
sition and fire control capabilities. These factors are normalised against a consistent 
standard or assigned a graduated evaluation score. The relative weights are based 
on the consensus judgment of experienced pilots as well as military analysts. The 
final performance are scores based on weighted calculations. The reader is referred 
to Hildebrandt, et al., 1986, and Congressional Budget Office, “Tactical Combat 
Forces of the USAF, Issues and Alternatives,” April 1985, for further details on the 
TASCOM methodology. 

The values of the parameters for future versions of each alternative are pro- 
jected as X;,q.,,. These are discounted (for deterioration) through the years of usage 
by the function d,(t). Technical expertise is required for such forecasting. Require- 
ment levels X;,,, need not be static but may be varied through the future years 
in accordance with projected threats. Figure 4.2 illustrates how this requirement 
level may change with time. Projections are required only for the span of the un- 
derlying network (Figure 3.1) and a steady state should be assumed beyond the 
forecast horizon. ‘These projections form the effectiveness aspect of DSRPM scenar- 


ios. Examples of U; { } and d;(_ ) are shown in Figures 4.3 and 4.4 respectively. 
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The utility function U;{ } is the instrument that allows us to express non-linear 
relationships such as that of diminishing marginal utility. Adjusting it regularly for 
the effect of changes in enemy force structure, size, tactics, strategy, etc., also serves 
to account for the interactive nature of defense “posturing”. The example shown 
here is of explicit utility functions and single mission systems. Implicit approaches 
may be used if utility functions can be specified explicitly only with considerable 
difficulty. Multiple mission systems may be modelled by combining utility values 


using weights or other multi-criteria modelling schemes. 
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Figure 4.2 c: Change in MOE Value and Requirement Level 


Figure 4.2: Illustration of Change in MOE Value and Requirement Level 
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Note : Utility functions need not be continuous. 
li MOEs for which low values are preferred to 
Utility . ; 
Value high values have decreasing functions. 





Value of MOE 


Figure 4.3 : Examples of Utility Functions 
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Figure 4.4: Examples of Deterioration Functions 


D. COST ESTIMATION 

Resource requirements are evaluated for both incumbent systems and prospec- 
tive alternatives. Cost models for the former deal with the projection of ownership 
costs based on actual cost experience as well as forecasts of disposal cost or salvage 
value. Those models developed for future alternatives are either statistical ones 


based on the physical and functional characteristics of these alternatives (Figure 


oO 


4.1) or analogy types based on comparison of similar systems and judgement con- 
cerning complexity factors. Procurement costs are included in cost estimation of 
alternatives but are considered sunk costs in the case of incumbent systems. 

Cost models used in DSRPM differ from the ones typically used in evaluations 
in that various types of costs such as procurement cost, operating cost, etc., are 
estimated not just for one possible life-cycle but for multiple possible ones in the 
future. The projection of costs needed to adapt cost models for use in DSRPM is 
usually accomplished using time-series or regression analysis. Adjustments may be 
made if significant and relevant technological breakthroughs are anticipated. For ex- 
ample, if major advances in certain automation technologies are expected five years 
from now, then manpower costs may be adjusted accordingly. Such adjustments 
should be made conservatively and in a not too abrupt manner. Figure 4.5 shows 
an example of the cost parameter forecasts that are used in the model. 

To examine how costs may be projected into the future, let us consider the 


following simplified life-cycle cost model: 
LC GME AC, aie je OC ee =) DC) (4.3) 
DAG MCs) — Ome, ae 2 ls (4.4) 

where: 


e LCC,,(t) is the life-cycle cost of alternative a if procured in year y and used 


for ¢ years 
e AC, is the acquisition cost of alternative a if procured in year y 
e MC,, is the maintenance cost of alternative a in the first year if procured in 


year y 


ol 
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Figure 4.5: Examples of Cost Parameter Forecast 


¢ OC, 18 the operating cost of alternative a in the first year if procured in year 


y 


e DC,,,(t) is the disposal cost of alternative a if procured in year y and used for 


t years (Note: This may be negative if asset has resale value.) 


Uay is the upper bound of DC,,,(t) 


Vay 18 a proportionality constant dictated by the resale market 


Ya,y> ay and ay are exponential constants 


O2 


If inflation and technological change have been (and presumably will be) gradual, 
the individual costs may be forecasted using non-linear regression models such as 
the following: 

ACuy = PAC,(1 + rac)” + ey (4.5) 


where: 


e PAC, is the present acquisition cost of alternative a 


@ rac is the compound inflation-technological improvement rate for acquisition 


cost 


e e, is the error term 


Solutions are obtained using either the maximum likelihood approach or Zellner’s 
method [Judge, et al., 1988] based on historical data. If major breakthroughs are 
expected in the production of alternative a, the forecasted figures for AC,,, may 
be adjusted to some extent. Anticipated foreign exchange rates may also be fac- 
tored into cost projections if they are relevant to a system renewal problem. The 
cost aspect of DSRPM is manifested in these projections. There are many differ- 
ent methods of forecasting both cost and effectiveness [Ayres, 1969, Stewart and 
Wyskida, 1987, and Ascher and Overhold, 1988], the choice of which depends very 


much on the degree of understanding one has of the dynamics of the change process. 


E. OTHER ELEMENTS 

The deterioration function d;(t) introduced in tlie section on effectiveness mod- 
elling is an example of how the change in performance parameters over the life of a 
system may be modelled. Audit programs such as OT&E, inspections and exercises 
as well as technical trials and laboratory tests provide the data needed to derive 


these functions. To see how deterioration may be incorporated in a cost model, let 
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us consider the following maintenance cost model: 
= = cht? (4.6) 
where: 

e M is the cumulative maintenance cost incurred 

e tis the age of the system 

e cis a proportionality constant 


e a and # are exponentiality constants 


This model links the rate of change in maintenance cost to the age t as well as past 
maintenance effort measured in terms of cumulative cost. Given an initial condition 
of M =0 at t = 0, we have: 

M =A (4.7) 


where: 


(4.8) 
A= <| a (4.9) 


al (4.10) 


Hence the rate of increase in maintenance cost is given by: 


— = Ayt?! (4.11) 


The parameters A and y are obtained by using regression analysis on historical 


maintenance cost date. A is the maintenance cost incurred in the first year of use. 


As 7 is always greater than 1, the rate of increase in maintenance cost ant will 
always be increasing over time. The impact of inflation and technological advances 
in maintenance cost is modelled by varying both A and 7¥ for future alternatives. A 
reduction in operating cost for future generations may also be modelled in a similar 
fashion. 

Operational requirements for present systems are projected by defining the 
minimum requirement level in each year for various performance parameters of a 
system. These are denoted as X ; ,,, in section B. For instance, shifts in enemy beach 
defence tactics may result in significant increases in the speed required of amphibious 
assault craft. Requirement projections form an integral part of a DSRPM scenario 
and are driven mainly by the net assessment process (Figure 4.1). 

Force levels decline even in peacetime with attrition due to accidents, wear 
and, in the case of human systems, retirement and resignation. If the number of an 
alternative a’ acquired in year y is taken to be the performance parameter Xj.’ y, 
the decline due to attrition can be modelled by the deterioration function d, ;(t). 
The effectiveness of various force levels can then be modelled using the associated 


utility function U; {}, and the force level requirement in the year (y +¢) is simply 
X 


X;y41- Treating the number acquired as a performance parameter in this way 


allows us to check the adequacy of the size of an acquisition. The attractiveness of 
a phased build-up versus the economies of a bulk buy may be studied by defining 
alternatives representing different systems and force levels. DSRPM _ produces an 
acquisition program that meets projected force level requirements by selecting and 
scheduling these alternatives. Figure 4.6 illustrates the kind of solution that may be 
generated by DSRPM in sucli studies. Force mix can be optimised by ensuring that 
the lagrangian multiplier associated with the levels of the various force coiponents 


are similar. In this way, force level and mix issues can be dealt with by the model. 
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Acqulsitlon In 3 lots — Initial lot of 6, 
then one of 12 for the bulld—up followed 
by one of 6 to moke up for attrition loss. 


Years 


Figure 4.6: Phased Acquisition Solution Example 


The last element we look at concerns the selection and scheduling of alterna- 
tives for implementation. The inputs for these are the cost and effectiveness models 
discussed in sections B and C of this chapter as well as budgetary, temporal, legal 
and other constraints. The combined dynamic programming — Lagrangian relax- 
ation discussed in Chapter 3 is applied to these inputs to produce solutions for the 
fixed budget and fixed effectiveness approaches. Constraints in the form of quantt- 
tative rules can be incorporated within DSRPM. An example is that of lead time 
requirements for various alternatives which is implemented in the prototype that 


was developed. The model is exercised by iterating it for various scenarios to make 
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explicit uncertainties faced by the decision maker. Strategic considerations with 
the formal quantitative results from DSRPM are then the basis of system renewal 
decisions which are the output of the force/resource planning process. 

DSRPM is easily implemented in phases. The model could first be applied to 
a couple of incumbent systems for which systems renewal decisions are expected in 
the next one or two years. This would allow enough time for a full implementation 
of the model for these systems. This phase serves to validate the model with a 
real application and the experience gained would guide implementation in the next 
phase. The model should however be applied as early as possible in a system’s 
life-cycle to provide enough time for the cost and effectiveness models to be refined 
with inputs from operational trials, exercises, wargames, etc. Therefore, the second 
phase should see the model applied to some new acquisitions. When management is 
comfortable with the model and is ready to introduce it on a broader scale, a master 


implementation plan can then be formulated to formalise the use of the model. 


Of 


V. PROTOTYPE DESIGN 


A. GENERAL 

To foster a better understanding of the DSRPM, a prototype was developed 
based on a hypothetical weapon system. This chapter documents the prototype 
which was developed on an IBM AT-compatible personal computer using PC SIM- 
SCRIPT II.5. The language was chosen for its superlative power and ease of use, 
as well as its provisions for incorporation of simulation in later development. The 
program is designed to produce for each run a single solution path corresponding to 
a point on Figure 3.2. The program reads data from an input file “IN.DAT” and 
channels results to an output file “OUT.DAT”. The type of solution depends on the 
approach selected for that run, namely Lagrangian relaxation, maximum effective- 
ness, minimum satisficing cost or maximum system cost effectiveness ratio. Solution 
approach is specified by the user in the input file. The budget is relevant only when 
the Lagrangian method is used. The overall structure of the prototype is described 


in Figure 5.1 and the program is listed in Appendix A. 


98 





RESULTS.PRINT 










CALCI.NODE 


--> denotes links 
conditional on 
input information 








PREAMBLE 













Print-out of 
intermediate and 
final results 


Declaration of 
entities and global 
variables 


Structuring the 
system renewal! 
network 












CALC2.NODIE 













Computing 
year & system 
ussociated with 

each node & 


ares entering it 


CALC.ALT.CODE 


INITIALIZE 





BEST.PATII 


Solution of 
optimal path 
using forward 

dynamic 
programing for 
given approach 































Read approach, 
scenario & other 
inputs. Create 
system renewal 
network & 
compute life-cycle 
effectiveness & 
cost of each arc 



















Initialize system 
renewal network 
& solve for optima! 
path given 
approach, scenario 
& other inputs 









Assigning 
alternative 
code to each arc 


ECHO.INPUT 


Check-print of 
input and key 
computational 
values when 
required 


















SUM.LCC.LCE. 


Compute 
cumulative life- 
cycle eff. & cost of 
3 systems in 
optimal path 














CALC.ARC 


Coinputation of 
life-cycle 
effectiveness & 
cost associated 
with each arc 




















Solves for optimal 
path & associated 
A* using 
supporting plane 
method 










Figure 5.1: Prototype DSRPM Block Diagram 
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B. ASSUMPTIONS 
The following assumptions are peculiar to the hypothetical weapon system 


upon which the prototype was developed: 


e The effectiveness sub-model is given by Equation 4.3 which was discussed as an 
example in the previous chapter. Deterioration functions d, ;(¢) are assumed 
to be exponential for all a and 7. Furthermore, utility functions U; {X} are 


assumed to be linear and bounded (maximum value of 1) as follows: 


xs 


Therefore, the total effectiveness of an alternative a acquired in year y in the 


t’® year of use is given by: 


TE, ye — > WU ardent (5.2) 
j=l 
= J iio _ posal! 


— Maxwve 
j=) a,y { J:@,y 


e The cost sub-model used is also similar to the example discussed in the previ- 


ous chapter. Annual maintenance cost AAfC,,;(t) is assumed to be: 


AM Cy y(t) = Kaytayt??! (5.3) 


At the end of a systein’s life, all parts that can be salvaged are sold and those 
parts that could not be sold are disposed of at a certain cost. Therefore, the 
model has both disposal cost and salvage value. Operating cost, disposal cost 


and salvage value are all assumed to vary exponentially through the years. 
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Therefore, the total life-cycle cost incurred by an alternative a acquired in 


year y and used for ¢ years ts: 


OC a (L) ar 1 Oe at (KOR Ce + OV ea 


+2 {Kay tayt"! + OC, ye°"'} (5.4) 
t=1 


C. MODULES 


This section briefly explains the function and workings of each component 
module of the prototype. 
1. PREAMBLE 
This module is peculiar to the SIMSCRIPT II.5 language. It provides 
for declaration of entities and global variables. Entities are structured data items 
that represent elements of a model. Definition of an element as a permanent entity 
entails the automatic creation of separate arrays for each declared attribute when the 
entity is created, the length of these arrays being the number of that element. For 
example, alternatives were defined to be the permanent entity ALT with attributes 
ALT.NAME and LEAD.TIME. When the entity ALT its created, two arrays with 
the names of the attributes will be created simultaneously and their length would be 
equal to the number of alternatives N.ALT (a system variable). Temporary entities 
would be relevant only in the modeling of dynamic systems where elements enter 
and leave. 
2. MAIN 
As the name suggests, this module controls tle overall computational 


flow, which is shown graphically in Figure 5.2. Note that the modules BEST.PATH 
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and SUM.LCC.LCE are executed several times if the Lagrangian relaxation ap- 
proach ts selected. 
3. INITIALIZE 
This module reads the input data, controls the modules generating the 
system renewal network and creates the permanent entities and associated arrays. It 
also reserves the memory space needed by the arrays. The flowchart for this module 
is shown in Figure 5.3. The input data is read in free format in the order shown in 
the listing. 
4. CALCI.NODE 
The structure of the network is determined in this module based on the 
number of alternatives and the minimum and maximum life-span of each system. 
The number of permanent entities NODE and ARC to be created is also computed. 
5. CALC2.NODE 
This module computes for each node information such as the year and 
system (transition) it belongs to, as well as the arcs entering it. The hierarchical list 
data structure illustrated in Figure 5.4 was used to represent the network. Therefore, 
the indices of the arcs entering node n runs from EP(n) to [EP(n +1) — 1]. 
6. CALC.ALT.CODE 
The code of the alternative represented by each arc is assigned in this 
module. 
7. CALC.ARC 
This module determines the tail node, life-cycle cost and effectiveness 
associated with each arc. The one-time cost items in Equation 5.4 are discounted and 
summed separately from the cumulative items such as maintenance and operating 


costs. Acquisition cost of the present system is ignored as sunk cost. The benefit 


accrued in each year is added to yield the cumulative or life-cycle effectiveness of 
the arc representing a given alternative and life-span. If the value of any MOE 
falls below that required in any year spanned by an arc, the life-cycle cost of that 
arc is set to an unfavorably high value to mark it as infeasible and to prevent its 
inclusion in the selection path. The highest index of the MOE(s) that failed to meet 
its requirement(s) is recorded as the limiting MOE in such cases. The flowchart for 


this module is shown in Figure 5.5. 
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START 


Call INITIALIZE 
Input control information 
Create network representing 
system renewal process 


Input projected cost and 
effectiveness parameters 
Compute life-cycle cost and 
effectiveness associated with 
each arc in the network 


Others Lagrangian 
Relaxation 
Call LAGRANGIAN 
@ Determine initial upper 
and lower bounds on A 
@ Determine A* by iterations 
of the supporting plane 
method 
@ Note: Fach iteration 
entails calling 
BEST.PATIHI and 
SUM.LCC.LCE to solve 
for L(A) and slope of LQ). 









Call BEST.PATH 
@ Uses dynamic program- 
ming to solve for path 
that is optimal under 
given approach 














Call SUM.LCC.LCE 
Compute cumulative 
life-cycle cost and 
effectiveness along 
optimal path 







Call RESULTS.PRINT 
Output of optimal solution for 
given approach based on 


given scenario (and budget in 
the case of the Lagrangian 
relaxation approach) 





END 


Figure 5-2: Flowchart or MAIN Program 
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Open input and output files 





Read input file 
Create permanent entities 
SYS, ALT & YEAR 
Reserve memory for arrays 
Call CALC1._ NODE to compute 
structure of system renewal network 


Create permanant entities 
NODE & ARC 


Assign alternative codes to 
representative arc 


Computes tail node, life-cycle 
cost & eff foreach arc 









Print title, pproach 
& budget 






CallECHO INPUT to 


print key inputs and No 
intermediate 
results 

Return 


Figure 5.3: Flowchart for INITIALIZE Module 
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TRANS. YEAR: ‘Transition year 


SYS.CODE: Preceding system code 
EP: entry potnt NODE 
TAIL: ‘T'ail node 
ALT.CODE: Alternative 
represented 


Note: The convention shown here is that 
of reverse star. 


Trans | Sys. 
Year | Code 
on] a 


= 
oO 
© 
O 


‘ 


pS 


Lo 
oe 
wo 
D> 


> 


Network Represented 


1993) 223.5 2 eee ee 
—— Alternative 0 


—— Alternative 1 


pe ya Alternative 2 





oO 
O 
NJ 


oO 
Oo 
WW 







a}ion| 
6 
oO 


1989 


Current 
Year System A System B 








Hierarchical 


List Data 
Structure 





ARC 


Figure 5.4: Hierarchical List Data Structure for System Renewal 


Network 


66 


Se 
CQ o 
o: 









START 





Initiahze life-cycle eff and 
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eff for various possible life- 


spans of the alternatives 
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Figure 5.5: Flowchart for CALC.ARC Module 
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Includes checks of 

MOE value (which 
deteriorates) ineach 
year to the level 
required in that year 
Failure to meet 
requirement for any 
MOE tn any year 
spanned by an arc results 
in the life-cycle cost 
being set toa very high 
value 


8. ECHO.INPUT 
Key inputs may be echoed for debugging purposes by setting the flag 
variable ECHO.ON in the input file to “on”. This module will then be activated 
and intermediate results such as tlie life-cycle cost and effectiveness associated with 
each arc are printed along with the key input data. 
9. LAGRANGIAN 
This module produces the optimal path for a given budget using the 
Lagrangian relaxation method. The value of A* corresponding to the optimal path 
is found by iterations of the supporting plane method. Each of these iterations 
entails calling the modules BEST.PATH and SUM.LCC.LCE to solve for L(A;) and 
the slope of L(A) at 4;. The limitations of binary representation will be exceeded, 
even with double precision variables, if 4; and/or the cost values are too high. The 
initial upper bound A,,4, given by Equation 3.18 works well for cost values scaled 
to no more than 10°. Figure 5.6 shows the flow of the computations in this module. 
Note that even though the slope at the minimum of L(\) itself may be positive or 
negative depending on whether A3 = A2 or \3 = A, respectively, the optimal solution 
path is given by the positively sloped line segment that is adjacent to the minimum 
point. When A3 = Aj, this line segment corresponds to A2 (otherwise, the intercept 
of the line segments at A; and 2 would not have yielded 43 = A,). Therefore, the 
slope is readjusted when A3 = Aj. 
10. BEST.PATH 
This module produces the optimal] solution path for the specified approach 
using forward dynamic programming. The criterion function L(J) associated with 
each arc J when the approach is to maximize the effectiveness to cost ratio of the 


individual systems is given by: 
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LCE(J) 


L)= (LOO@)= MIN LCC 41) 


(5.5) 


where: 


e LCC(J) and LCE(J) are the life-cycle cost and effectiveness respectively of 


are J 


e MIN.LCC = min LCC(J) 


This is to allow for the possibility of LCC(J) being zero or negative. 
Provisions are also made to ensure that the solution path has sufficient lead time 


for each alternative chosen. The computational flow is illustrated in Figure 5.7. 
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Figure 5-6: Flowchart for LAGRANGIAN Module 
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Figure 5.7: Flowchart for BEST.PATH Module 
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11. SUM.LCC.LCE 
The solution path is traced backward from the terminal node with the 

highest value of the criterior function CRIT.FUNC(]) and the life-cycle cost and 
effectiveness of the systems are summed along the way in this module. Note that 
all the terminal nodes will have very low values (negative) of CRIT.FUNC(I) if all 
paths violate either lead time or MOE requirements or both, i.e., no feasible solution 
exists. Conversely, if CRIT.FUNC(I) of a terminal node is not a very low value, 
below 107°, then a feasible solution path exists that leads to that node. 

12. RESULTS.PRINT 


Computational results are printed by this module. 


D. COMPLEXITY AND EXECUTION TIME 

In terms of time complexity, the dominating piece of code is the CALC.ARC 
module. Let the difference between the maximum and minimum life-span of the 
present system be n, and that of successor systems be n2,73,---,3, where s is the 
number of stages or systems in the systems renewal network. If we assume n — n; 
for Vi > 1, then the structure of the network is given by ‘Table 5.1. 


The total number of arcs n is thus given by: 


s—2 
n = ltn,+a(s—1)(1+m)(ltn)+an(l+n)>>i 


t=] 


l+n,+a(s—1)14+n,)1+4+n)+an(1+n)(s —2)(s—1)/2 (5.6) 


If we let the average life of successor systems be denoted as 1, and the number of 
MOFEs evaluated be e, then the growth rate of execution time, ET, 1s of the order 


m(l +e). In other words, 


TABLE 5.1: System Renewal Network Structure (n,; =n for V; > 1) 
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And so on until 2 = s 











Note: a denotes the number of alternatives evaluated 


ET = 0O|m(l+e)] 


0 {1 +n t+a(s—1)1+n,)(1+n)+ eee cr) 


0 |max(an, ns,an*s*)max(l, e)| (oe) 


Therefore, execution time for the prototype DSRPM is of polynomial order, the key 
determinants of which are n and s. We can also see that time complexity is linear 
iigcwmeanGee.ror a run with @ = 2, = 5o5n.—.5,5 = 3.1 = 7.5.and e =3 "the 
execution time of the prototype DSRPM on the Zenith AT-compatible ts typically 


about 43 seconds. 


VI. CASE STUDIES 


Case studies were conducted using the prototype to demonstrate the plausi- 
bility and usefulness of the DSRPM. To better serve the purposes of illustration and 
plausibility testing, the case studies were designed to comprise a series of scenarios 
of incremental complexity, all of which are based on the same hypothetical system 
discussed in Chapter V. Although effort was made to incorporate as much realism as 
possible, the case studies and the weapon system were kept simple to better portray 
the workings of the model. A good understanding of and confidence in the model is 
necessary before one confronts the complexities involved in a real application. 

In the case studies, the key effectiveness parameters or MOEs are assumed to 
be range (nautical miles), accuracy (percentage hit) and availability (percentage) 
with relative weight IV; of 0.4, 0.3 and 0.3 respectively. The cost and effectiveness 
characteristics of the incumbent system assumed are shown in ‘Table 6.1. Disposal 
cost is assumed to be zero for the incumbent system and alternatives, i.e., DC, = 0 


fOr Von 


TABLE 6.1: Characteristics of the Incumbent System 
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As astarting point, we shall assume that two alternatives, 1 and 2, are available 
and that their minimum and maximum life are five years and ten years respectively. 
They are assumed to require no development or acquisition lead time (presumably 
available “off-the-shelf”). The incumbent system (alternative 0) is assumed to have a 
minimum and maximum life of zero years and five years respectively, i.e., alternatives 
may be introduced immediately, or at the latest, five years hence. 

Four solution approaches are used in each of the case studies, namely the 
minimum average annual cost (MAAC) approach, the minimum annualised cost 
(MAC) approach, the system effectiveness to cost ratio (SECR) approach and thie 
Lagrangian relaxation (LR) approach. The average annual cost is obtained by di- 
viding the total cost of a succession of three systems (in current dollars) by the years 
spanned by the systems. The MAAC approach thus ignores monetary interest. The 
MAC solution is obtained by amortizing the total cost of the system over the years 
spanned. It is the solution obtained when one wishes to minimise cost needed to 
meet effectiveness requirements. The SECR approach maximizes the cumulative 
total of the effectiveness-to-cost ratio of each of the three systems in a path. The 
LR approach was explained in Chapter 3. Solution points obtained using the LR 
approach are linked as in Figure 3.5 to form the efficiency boundary. The point on 
the boundary that yields the maximum overall effectiveness to cost ratio would be 
the preferred solution if one adopts the approach of maximising returns on cost. 

A notational point to note is that the transition years shown in the results of 
the case studies are all abbreviated to the last two digits. That is, 97 and 05 mean 


1997 and 2005 respectively. The current year is 1989. 


—T] 


en 


A. CASE STUDY 1: STATIONARY ENVIRONMENT 

The first scenario we examine is one in which the only change occurring is 
that due to deterioration. Technology and requirement levels are at a standstill and 
no inflation occurs. The value of money is constant over time, 1.e., interest rate is 
zero. This highly improbable scenario is implied by the commonly used “equal life” 
assumption whereby an alternative, once selected, will be selected over and over 
again and will always have the same service life. 

Supposing alternatives 1 and 2 have the following cost and effectiveness char- 
acteristics shown in Table 6.2 and 6.3. Alternative 2 is thus three times as expensive 
to acquire but has a longer range and offers lower maintenance and operating costs 
than Alternative 1. Requirement levels are assumed constant at 95 nm for range, 
80% for accuracy and 70% for availability. Results of analysis performed using the 


prototype DSRPM for this case is presented in Table 6.4 below. 


TABLE 6.2: Characteristics of Alternative 1 






[MOE [Vahie When New[ @ [Cost Parameters. —SSCi*d 
(Range(am) | 120 | -0.03 | Acquisition [Maintenance | Operating | Salvage | 
[Recuracy (%) [95 —C—~O| AG Th] Joc, & [Sh] af 
“Availability (| __95 [0.03 100 [30.1.3 | [0.05] 0/02) 






The value of the total life-cycle cost and effectiveness are plotted in Figure 6.1. 
As interest rate is zero, MAAC and MAC approaches yielded the same solution. The 
results show that an alternative that is selected will be selected repeatedly in the 
future, the service life each time being the same. This is only to be expected as, 


under static conditions, there are no reasons for change. Note that the solution that 
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TABLE 6.3: Characteristics of Alternative 2 












[MOE [Value When New| [~~ Parameters id 
[Range (am) | _140 | -0.03| Acquisition [ Maintenance | Operating | Salvage_| 
Accuracy (%) | 98_—~+4f-002) AC? [el [OCG] & [Sta] 0 | 
[Availability (%| 95 | -003) 300 fi5[ 13] 4 [0.05 -90 [0.2] 





TABLE 6.4: Results for Stationary Environment 








[ Solution | Total LOC | Total LCE [LCE/LCC| Att. Chosen | Transition Vear| 7] 
[Approach |(3 Systems) | (3 Systems)[ Ratio [1[2] 3 [1]2] 3 | | 
MAAC | 120605] 16580] 0.0i372|0]2| 2 |9ijoo| 07 |_| 
[MAC [1206.05 16.550] 0.01372 [0/2] 2 |o1[99] 07 |_| 
[SECR [1097.35] _15.002[ 0.01373 /0/2| 2 [s0/o7| os |_| 
[in [1270.81] 17.265 [0.01359| 0/2| 2 [80/00 08 | 0 _| 
[IR | 1206.05 | _16.550| 0.01372 /0/2| 2 [91}99| 07 | 00110 | 
tr | _inar.33] 15.816] 0.013790 [2] _2 [90/98] 0G | 0.01251 | 
(ER | 801.8i] 10.934 | 0.01367 0]1 [1 [90] 95] 00 _| 0.01413 
PIR | _731.86] 10179 | _001354) 0), 1 [89/94] 99 [0.01509 











gives the maximum return per unit cost for the three systems overall (marked with 
an asterisk in Table 6.4 and a square in Figure 6.1) is different from that obtained 
using the SECR approach. The overall LCE to OCC ration in the case of the SECR 
solution is not as high even though the solution also lies on the efficiency boundary. 
The implication of this is obvious —- maximizing return per unit cost for individual 
systems considered separately leads to a good but not best solution in terms of 


overall return per unit cost. 


—~] 
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Efflclency Boundary 
(formed by Logranglan Solution Points) 
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Figure 6.1: Stationary Environment Case Study Results 


B. CASE STUDY 2: INFLATION AND INTEREST RATES 


The next scenario we look at is probably the one most commonly encountered 
in a typical system renewal study. Future inflation and interest rates are assumed 
to be non-zero and constant, and future technological advancements and changes 
in requirement levels are ignored. For illustrative purposes, let us say the inflation 
and interest rates were projected to be constant at four and five percent per annum 
respectively. Cost projections based on these rates are shown in Figure 6.2 in the 
sample input file for this case study (Appendix B). Analysis based on this scenario 


produced the results shown in Table 6.5. 
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Figure 6 2b: Cost Evolution for New Alternative 2 


Figure 6.2: Evolution of Cost Parameters for Alternatives 1 and 2 
With Inflation (No ‘Technological Progress ) 
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TABLE 6.5: Results for Inflation and Interest Rates 


[ Solution | Total LCC | Total LCE LCE/LCC [ Alt. Chosen | Transition Vear| 4} | 
[Approach [8 Systems [3 Systems |" Rafi] 3 [ata] 3 | 


1031.26 0.01674 
§11.98 0.01647 
921.03 0.01563 
1031.26 ip 0.01674 
986.44 ; *0.01678 
944.29 0.01675 
§11.98 Oto 0.01647 
683.34 0.01600 
643.73 sii 0.01581 


0 
0.01596 
0.01742 
0.01845 
0.01898 
0.01904 





mm Fr KD RD BHD HO LO LY 
mm mm Mm KD AD RD BR 


0 
0 
0 
0 
0 
0 
0 


Total LCC and LCE from this table are plotted in Figure 6.3. The efficiency 
boundary is shifted towards lower cost since the interest rate is higher than the 
inflation rate. An opposite movement would occur if the converse was true. Note 
that the shift is greater for the higher cost end because inflation and interest rates 
have a bigger impact there. As the interest rate is significantly non-zero, we see 
that the solutions obtained using the AfAAC and Af AC approaches are different. 
Note also that the SECR solution is clearly inferior as it lies somewhat below the 
efficiency boundary. From the above table, we can see that Alternative 1 is the 
low-end alternative. The last line in the table shows the solution that meets the 
specified minimum requirements at the lowest possible total cost, which is the one 


whiere the service lives of all three systems are at their minimums. 
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Figure 6.3: Inflation and Interest Rate Case Study Results 


C. CASE STUDY 3: TECHNOLOGICAL ADVANCEMENT 

For this scenario, let us suppose that technological progress has the effect of 
suppressing cost growth due to inflation from four to two percent per annum. The 
revised cost projections are shown in Figure 6.4. The sample input file for this case 
study is shown in Appendix B. In addition, technological advances are also reflected 
in improvement in performance characteristics. For the purpose of this case study, 
let us suppose that the various MOEs evolve as shown in Figure 6.5 and the sample 
input file in Appendix B. Improvements are primarily in range, with Alternative 1 
becoming marginally better than Alternative 2 about ten years from now. Other 


conditions are as assumed in Case Study 2. The results for this case study are shown 


in Figure 6.6 and Table 6.6. 


Si 


Contrasting ‘Table 6.5 and Table 6.6, we note that Alternative 1 is now the 
dominantly favored alternative because it became competitive performance-wise as 
well as cost-wise. Technological advancement resulted not only in a further shift 
in the efficiency boundary towards lower cost but also an extension into higher 
levels of cumulative effectiveness. The lower end of the boundary is not affected 
as much as solutions in the higher end that extends further into the future. Note 
that technological advancement also results in lengthening of the maximum service 
life (three systems), as evident in the maximum effectiveness solution identified by 
A* = 0. In Case Studies 1 and 2, no solutions extend beyond the year 2008 whereas 


in this case study, some solutions extend to the year 2009. 


TABLE 6.6: Results for Technological Advancement 


[Sill [Total LOG Wel | ESE /LCOT NTE ges gree eal | 
[Approach | Systems [3 Systens |" Racio [if2] 3 ji [2 {3 | 


0.02082 
0.01971 
0.02015 


0.02068 
*0.02082 
0.01926 
0.01889 
0.01828 


0) 
0.01444 
0.02473 
0.02601 
0.03180 


eet OAD Or 
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Figure 6. 4b: Cost Evolution for New Alternative 2 


Figure 6.4: Evolution of Cost Parameters for Alternatives l and 2 
With Inflation and ‘Technological Progress 
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Note : Technological progress is primoarltly In 
the 1st MOE for both citernatives. Improvement 
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Figure 6.5: Evolution of MOEs With Technological Progress 
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Figure 6.6: Technological Advancement Case Study Results 
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D. CASE STUDY 4: REQUIREMENT ESCALATION 
Let us suppose that everything is as before, except that requirement levels set 
for range and availability increase over time as shown in Figure 6.7. The results for 


this scenario are shown in Table 6.7 and illustrated in Figure 6.8. 


TABLE 6.7: Results for Requirement Escalation 


[Seti [Tod LOC [Tova LOE [LCELCC AN, Chosen Traction Ver] OE] 
Approach | 8 Systems| Systems |" Ratio [fa] 3 [aya y 3 | 


0.01955 
0.01955 
0.01926 


0.01883 
+*0.01961 
0.01955 
0.01926 
0.01828 


0 
0.00097 
0.02011 
(02729 
0.02863 


— | S| SS RDS eS 





Comparing Figures 6.6 and 6.8, we can see that the escalation in requirement 
levels resulted in a downward depression of the higher cost end of the efficiency 
boundary. This happens because the escalation in requirement levels results in the 
shortening of service lives only if the solution path extends far enough to be affected. 
Note that the maximum effectiveness solution, for instance, now extends only until 
the year 2006 whereas in the previous case study it reaches the year 2009. An 
interesting result in the equality of the AfAAC and AfAC solutions even though 
the interest rate is five percent per annum. Note also that solutions calling for the 


first transition to be made in 1991 are clearly dominating. 
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Figure 6.8: Requirement Escalation Case Study Results 
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BE. CASE STUDY 5: SCHEDULING OF UPGRADES 

Supposing Alternative 1 is an upgrading option, the exercise of which is sig- 
nificantly cheaper if carried out along with major scheduled overhauls. If suc: 
maintenance is carried out every three years, (e.g., refit of a ship), and the savings 
is say twenty percent, then the acquisition cost of a new Alternative 1 varices as 
illustrated in Figure 6.9. In this scenario, everything else is as assumed for Case 
Study 4. 

The results for this scenario are given in Table 6.8 and Figure 6.10. As evident 
from the table, there are no saving realised in the higher cost end where the discounts 
are not taken up. The lower end of the efficiency boundary shifts to the left as a 
result of the savings. We sce that the solution offering maximum effectiveness per 
unit cost suggests that the incumbent should be replaced with Alternative 2 in 
1990, then upgraded with Alternative 1 in 1997 (a non-maintenance year with no 
discount) for use until the year 2005. The SECR solution also indicates an upgrade 
carried out in a non-maintenace year (1998). The overall desirability of matching a 
major upgrading program to some maintenance schedule to miniinise downtiine or 


to secure some savings can thus be evaluated using the model. 
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Figure 6.9: Evolution of Cost Parameters With Scheduled Overhauls 
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Figure 6.10: Scheduling of Upgrading Case Study Results 
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TABLE 6.8: Results for Upgrade Scheduling 


[ Solution [Total LOC | Total LCE [LCE/LCC| Alt. Chosen | Transition Vear| Xj] 
[Approach | § Systems [3 Systems [Ratio [ia] 3 [fay 3 [| 


0.01955 
0.01949 
0.01926 


0.01883 

0.01961 0.00097 
*0.01962 0.01953 
0.01953 0.01981 





Note: In this case study, a 1 for alternative chosen represents an upgrade of the 
preceding system. Thus the sequence 0,2,1 means the incumbent is replaced with 
alternative 2 and then upgraded with alternative 1. 


F. CASE STUDY 6: DEVELOPMENT LEAD TIME 

Let us suppose that everything is as in Case Study 5 except that Alternative 
1 is an upgrade that requires one year of development and cannot be introduced 
immediately. The impact of such lead-time requirements can be seen in Table 6.9 
and Figure 6.11. From the table, we can see that the lowest cost solution (bottom of 
table) still involves an immediate transition as in the other case studies. However, 
the transition is to Alternative 2 in this case as the model now bars the immediate 
execution of Alternative 1. The diagram clearly shows that the lead-time require- 
ment pushes up cost at the lower cost end of the efficiency boundary. The MAAC 
solution, the SECR solution and the maximum effectiveness per unit cost solution 
are all unaffected. The MAC solution, however, became equal to the MAAC solu- 


tion (curiously). Note that the MAC solution in Case Study 5 involved immediate 
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introduction of Alternative 1 (upgrade of alternative 0 preceding it); a solution that 


is not admissible in this case study. 


TABLE 6.9: Results for Development Lead-Time 


[Salta [Total LOC Totel LCE TLE/LGG | A Chosen [Tranation Veer] 
Approach | 3 Systems [3 Systems|" Natio [1[aT 3 [i fay 3 | 


a 01955 
0.01955 


0.1926 
0.01883 


0 
0.00097 
0.01953 
0.02063 
0.02729 
0.11078 


0.01961 
*0.01962 
0.01907 
0.01926 
0.01731 
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Figure 6.11: Development Lead Time Case Study Results 


G. GENERAL OBSERVATIONS 

From the series of six case studies analysed using the prototype DSRPM, we 
saw the impact of various issues as they were integrated into the model. DSRPM of- 
fers a common framework for decision makers to explore the consequences of various 
possibilities. It does not advocate any single solution for a given scenario, but rather 
maps out the efficiency boundary spanning the entire solution space and indicates 
the solutions obtained using various approaches. A plot of the results for various 
scenarios, as is done in Figure 6.12 for the case studies, gives the decision maker a 
good overview of the problem. The picture can be enriched by varying the shade of 
the result for each scenario according to its likelihood of occurrence. Furthermore, 


a third axis indicating overall service time for the three systems can be added as 


wi 


in Figure 6.13 to give fuller expression to the problem. As explained in Chapter 2, 


the system renewal problem is basically a three-dimensional one - cost, effectiveness 







and time. 
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Figure 6.12: Consolidated Results for Case Studies 


Although three systems are analysed in the model, the focus is really on when 
the first transition should take place and what alternative should be selected to 
succeed the incuinbent system. Analysis may reveal a dominating (good under 
all scenarios) or robust (favored under most scenarios) solution to the timing and 
choice involved in the first transition. For example, the escalation of requirement in 


Case Study 4 resulted in hastening of the first transition by a year for the maximum 
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Figure 6.13: 3-D Plot of Efficiency Boundaries 


overall effectiveness to cost ratio solution. In Case Study 6, however, we saw that this 
solution was not affected by the development lead time of one year. By varying the 
scenarios for various requirement escalation trends, development lead time estimates, 
etc., we can test the sensitivity of various solutions to variations in estimates and 
trends. Even if analysis proved inconclusive, a rough idea of how long the incumbent 
is likely to be kept will emerge. No commitment needs to be made if the time to 
the likely date of the first transition is more than the lead time required to effect 
the likely transition. For example, if analysis shows that an EW suite upgrade of a 


frigate should take place three years down the road and such an upgrade has a lead 
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time of two years, we needn’t make a commitment yet, but we know the decision 
is due a year from now. At that time, the model can be run again with updated 
information. In this way, we can ensure decisions are made in a timely, fashion using 
the latest information. 

Besides facilitating timely, consistent and rational decision-making, the model 
also looks forward to the future to provide estimates of cost, effectiveness and service 
life that facilitates good planning. If aggregated to sufficiently high levels, the model 
allows future year procurement budgets to be forecasted. Designs and support plans 
for succeeding systems can be optimised to service lives that are more realistic 
in terms of broad economic, technological and operational outlook. Amortization 
periods would no longer need to be arbitrarily assumed and budgets can be better 
planned knowing cost estimates of succeeding systems. The anticipated system 
selections also serves to better guide planning for defense infrastructure, R&D and 
other strategic concerns. The explicitness of the model fosters better consensus 
building and coordination among the various agencies responsible for system renewal 
planning. As discussed in Chapter 4, Section E, force level and mix issues as well 
as phased build-up versus bulk procurement decisions can also be dealt with by 
the DSRPM. The model thus provides a comprehensive and systematic framework 
for system renewal decisions to be made optimally and in a consistent and timely 
manner based on current forecasts. 

Lastly, it should be noted that inflation and interest rates were assumed con- 
stant and linear utility functions were used in the case studies for reason of conve- 
nience and not necessity. ‘The model allows these rates and effectiveness requirement 
levels to be varied for future years in a scenario. Long-range projections can be 
varied more than near-term ones to account for greater uncertainty in these projec- 


tions. The impact of technological advancement is also freely interpreted in terms 
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of changes in cost and effectiveness parameter projections and no restrictions are 
placed on the definition of utility functions. The DSRPM thus affords very flexible 


means to model the factors influencing systems renewal problems. 
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VIL. CONCLUSIONS AND 
RECOMMENDATIONS 


A. CONCLUSIONS 

The dynamic system renewal planning model (DSRPM) provides a system- 
atic framework within which optimal system renewal decisions may be made in a 
consistent and timely manner, under realistic conditions of changing environments. 
Conceptually, it embodies a concept of obsolescence that is more comprehensive, 
better formalised and more relevant to military planners than is available in current 
literature. Prediction of what future renewals are likely to be made and when, and 
estimation of the associated cost and effectiveness would allow plans to be made 
anticipatively with a long-term perspective. Through simple reformulation and use 
of shadow prices, bulk versus phased procurement decisions as well as force level 
and mix issues can also be examined using the model. 

The plausibility and usefulness of the model was demonstrated by a scries of 
case studies using the prototype that was developed. The equal life assumption 
often used in system renewal models was illustrated by the first case study. Other 
case studies showed the imnpact of interest and inflation rates, technological progress, 
requirement escalation and development lead time. Application of the model to the 
scheduling of upgrades that are cheaper to carry out during overhaul periods was 
also shown. 

The prototype DSRPM illustrates how the model may be used as a descriptive 
decision aid to better understand the complex effects of a diverse range of issues on 
the system renewal decision. It was shown how the model organises financial, oper- 


ational, intelligence and other staff projections within an integrated framework to 
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provide decision makers with a broad, consistent, long-term perspective of relevant 
issues. The decision maker is presented with a graphic picture of the entire solu- 
tion space as structured by the system renewal scenarios considered. Scenarios are 
flexibly modelled and are used to point to solutions that are robust in the face of un- 
certainty. By continually refining scenarios based on lessons learned from exercises, 
audits, and other new information on trends, past experience is brought to bear on 
the present so that the basis of decisions is kept current and realistic. The commu- 
nication among all parties involved in this process is facilitated by the explicitness 
of the model. This in turn fosters consensus building necessary for cohesive effort 
in defense planning. 

The model that was developed is a conceptual prototype that is as yet untested 
by a real application. Much remains to be done before it is ready for Be aiicatibii: 
wide implementation. However, it holds much promise and should be examined 


further for future implementation. 


B. RECOMMENDATIONS 
The prototype that was developed is a rudimentary one that could be extended 
in many ways as future research topics. It is recommended that research proceed in 


the following order of tractability: 


1. Factoring foreign exchange rates into the model is a desirable goal. The treat- 


ment of interest and inflation rates would, however, require some adjustments. 


2. Optimal force mix can be determined if shadow prices of different components 
of a task force are made compatible. Means of integrating the utility functions 


of these components under a common mission structure needs to be examined. 


. Simulation and forecasting capabilities can be incorporated to better explore 


the range of system renewal scenarios that may be encountered in the future. 


. Stochastic models can be constructed to better portray uncertainty in cost, 


effectiveness and requirement levels. 


. Utility modelling may be facilitated by the use of implicit approaches that do 
not require utility functions to be explicitly defined over the entire domain of 


the MOE value. 


. System renewal models developed for various systems could be linked and ag- 
gregated to higher entities using schemes that incorporate organisation struc- 
ture into the solution algorithm to avoid sub-optimisation. Plans developed 
for ships could, for example, be tied to plans for their operational task forces 
and administrative units. The ideal would be to achieve a level of aggregation 
high enough to interface with arms race and econometric models so that very 


high levels of optimisation may be achieved. 


. Fuzzy set and logic theory may be introduced to better model elements that 
are inherently vague. Requirement levels, for example, may be made elastic by 
using What is known as fuzzy numbers. This may allow some form of abstract 


qualitative reasoning to be represented in the model. 


APPENDIX A: PROGRAM LISTINGS 


1 preamble 

2 

3 normally, mode is real, dimension is 0 

4 

5 permanent entities 

6 

7 every YEAR has a DISCOUNT.FACTOR and a CUM.DF 

8 

9 normally, mode is integer 

10 

11 every NODE has a TRANS.YEAR, a SYS.CODE, 

i2 a BEST.PRED.NODE and a BEST.PRED.ARC 

13 every ARC has a TAIL and an ALT.CODE 

14 every SYS has a MIN.LIFE, a MAX.LIFE, a NUM.NODES, a NUM.ARCS and 
15 a BASE.YEAR 

16 every ALT has an ALT.NAME and a LEAD.TIME 

7 

18 define MOE, MOE.DF as 3-dimensional real arrays 

19 derinesne, MC.P1, MC.P2, OC, OC.EF, DC, DC.EF, SV, SV.EF and 

20 MIN.MOE.REQ as 2-dimensional real arrays 

Zell define MOE.NAME as a 1-dimensional text array 

22 define L and CRIT.FUNC as 1-dimensional double arrays 

23 define LCE, LCC, MOE.WT, O.MOE and O.MOE.DF as 1-dimensional 

24 real arrays 

25 define MAX.PRED.NODES, LIM.MOE and EP as 1-dimensional integer arrays 
26 define ARGMAX.CF, ARG.MAC, ARG.AAC, CURR.YEAR, NUM.MOE and NUM.ALT as 
27 integer variables 
28 define LAMDA, LAMDA1, LAMDA2, LAMDA3, SUM.LCC, SUM.LCC.1, SUM.LCC.2, 
7%) SUM.LCE, SUM.LCE.1, SUM.LCE.2, MAX.LCE, MIN.LCC, MOE.VALUE, 
30 SLOPE1, SLOPE2, SLOPE3, Ll, L2, L3, BUDGET, AMC, AOC, DCY, 
31 BVY OMe) Pl, O.MCoP2, ©.0C, 0.0C.EF, 0.DC, O.DC.EF, 0.SV, 

32 O.SV.EF, MAX.MAC, MAX.AAC, AAC.CF, INT.RATE and INT.DF as 

g) 8) real variables 

34 define MAX.CF as a double variable '' Double means double precision. 
35 define APPROACH, ECHO.ON, CASE.NAME and ALT.NAME as text variables 
36 

37 end 


she) 


main 


call INITIALIZE 


if APPROACH = "Lagrangian" 
call LAGRANGIAN 
else 


call BEST.PATH 
Call SUM. Dee- LCE 
10 always 


ODN AUS WN He 


eZ call RESULTS.PRINT 


14 stop 


100 


routine INITIALIZE 


open 1 for input, name is "IN.DAT" 
use 1 for input 

open 2 for output, name is "OUT.DAT" 
use 2 for output 


read CASE.NAME 


read N.SYS, N.ALT, NUM.MOE, CURR.YEAR, INT.RATE 
read ECHO.ON 


create every SYS and ALT 


N.YEAR = 1 
for I = 1 to N.SYS 
do 
read MIN.LIFE(I) and MAX.LIFE(I) '' Computing number of years to 
N.YEAR = N.YEAR + MAX. LIFE(I) '' be covered in the analysis. 
loop 


create every YEAR 


for I = 1 to N.ALT 
read ALT.NAME(I) and LEAD.TIME(I) 
read DISCOUNT. FACTOR(1) 
CUM.DF(1) = DISCOUNT. FACTOR(1) 
for I = 2 to N.YEAR 


do 
read DISCOUNT. FACTOR (I) 
CUM.DF(I) = CUM.DF(I-1) * DISCOUNT. FACTOR(I) 
loop 
read APPROACH '' Determines the type of solution to be generated 
read BUDGET '' Used only in the Lagrangian approach 


'' Allocating memory to arrays 

reserve MAX.PRED.NODES(*) as N.SYS 

Beserve AC(*,*), MC.P1(*,*), MC.P2(*,*), OC(*,*), GCC.EF(*, “3 
DC(*,*), DC.EF(*,*), SV(*,*), SV.EF(*,*) as N.ALT by N.YEAR 

reserve MOE(*,*,*) and MOE.DF(*,*,*) as N.ALT by NUM.MOE by N.YEAR 

reserve MOE.NAME(*), MOE.WT(*), O.MOE(*), O.MOE.DF(*) as NUM.MOE 

reserve MIN.MOE.REQ(*,*) as NUM.MOE by N.YEAR 


call CALC1.NODE '' Structuring the network. 
create every NODE and ARC 

reserve CRIT.FUNC(*) as N.NODE 

reserve EP(*) as (N.NODE + 1) 

reserve LCE(*), LCC(*), LIM.MOE(*) and L(*) as N.ARC 


call CALC2.NODE '' Computing the year & system of each node 
'' &€ the start index of arcs entering it. 


call CALC.ALT.CODE '‘'' Assigning code of alternatives. 


read O.MC.Pl, O.MC.P2, 0.0C, 0.OC.EF, 0.DC, O.DC.EF, 0.SV, O.SV.EF 


10] 


for I = 1 to N.ALT 
for J = 1 to N.YEAR 
read AC(1I,J), MC.P1I(1I,d), MC2P2(1,3),)OC(1,o)mmOC-Fr( 1.1) 
pc(I,J), DC.EF(I,J), SV(1I,J) and SV.EF(I,J) 
for J = 1 to NUM.MOE 
read MOE.NAME(J), MOE.WT(J), O.MOE(J) and O.MOE.DF(J) 
for J = 1 to NUM.MOE 
for K = 1 to N.YEAR 
read MIN.MOE.REQ(J,K) 
for I = 1 to N.ALT 
for J = 1 to NUM.MOE 
for K = 1 to N.YEAR 
read MOE(I,J,K) and MOE.DF(I,J,K) 


call CALC.ARC '' Computes tail node, LCC and LCE for each arc. 


print 3 line with CASE.NAME thus 
Dynamic End-Of-Life-Cycle Planning Model 


Case OCF KHKKHKKKEKKKEKKKKEKKKKAKAKAKKKKKAKKAKKKKKKE 


skip 1 line 
print 4 line with APPROACH and BUDGET thus 
Approach — #k&kkkkkkrrkkakkake Budget = $kkkkkakkn ke 


Note : Budget is used only for the Lagrangian relaxation approach and is 
ignored in the other approaches. 


if ( ECHO.ON = "On" ) 
call ECHO.INPUT 


skip 1 line 
print 9 line thus 
Candidate Arcs For Solution Path 


Note - LCC is made to be unfavorably high if the minimum requirement 
level for any component MOE of LCE has not been met. The limiting MOE 
is given as the highest index of the MOE(s) that failed to meet the 

minimum requirement level(s) during the life-span covered by the arc. 


Arc Alt. Tail Head Life-Cycle ' Limiting 
No. Code Node Node Effectiveness Cost MOE 
skip 1 line 
for I = 2 to N.NODE 
for K = EP(I) to ( EP(I+l) - 1 
print i line with K,UALT.COBE(K), TAIL(K), 1, LCE(K), LCC(rm 
and LIM.MOE(K) thus 
kak kak kak kik kahkhakh  kaak kahkhh kh kak 
always 
return 
end 
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1 routine CALC1.NODE 

2 

3 '' Structuring the network by determining the number of 
4 ‘' nodes and arcs in each system 

5 

6 BASE.YEAR(1) = CURR.YEAR + MIN.LIFE(1) 

7 MAX.PRED.NODES(1) = 1 

8 NUM.NODES(1) = MAX.LIFE(1) - MIN.LIFE(1) + 1 

9 NUM.ARCS(1) = NUM.NODES(1) 
10 N.NODE = NUM.NODES(1) + 1 

nel N.ARC = NUM.ARCS(1) 

le for I = 2 to N.SYS 

13 do 

14 I= 1 — 1 

15 BASE.YEAR(I) = BASE.YEAR(J) + MIN.LIFE(I) 

16 MAX. PRED.NODES(I) = NUM.NODES (I-1) 

i NUM.NODES(I) = NUM.NODES(J) + MAX.LIFE(I) - MIN. LIFE(I) 
18 NUM.ARCS(I) = N.ALT * NUM.NODES(J) * 

19 ( MAX.LIFE(I) - MIN.LIFE(I) + 2 ) 
20 N.NODE = N.NODE + NUM.NODES(T) 
21 N.ARC = N.ARC + NUM.ARCS(T) 

22 loop 
23 
24 return 
25 
26 end 
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rou 
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4 
> 
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10 
11 
12 


35 end 


tine CALC2.NODE 


‘* Computes for each node information such as the year and system 
'* it belongs to, as well as the lead index ( contained in entry 
‘' point array EP ) of the arcs entering the node 


TRANS.YEAR(1) = CURR.YEAR 
SYS.CODE(1) = 0 

EP(1) = 0 

K = 2 

JJ = 1 

for J = 1 to NUM.NODES (1) 
do 


TRANS.YEAR(K) = BASE.YEAR(1) + J - 1 
SYS.CODE(K) = 1 
EP(K) = EP(K-1) + 1 
K = K + 1 
loop 
for I = 2 to N.SYS 
do 
for J = 1 to NUM.NODES (TI) 
do 
TRANS. YEAR(K) 
SYS.CODE(K) = I 
EP(K)) = EP(K-1) + JJ 
M = MAX.LIFE(I) - MIN.LIFE(I) + 1 
JJ = N.ALT*min.f(J,MAX.PRED.NODES (I) ,M, (NUM.NODES (I) -J+1)) 
K = K +1 
loop 
loop 
EP(N.NODE+1) = EP(N.NODE) + N.ALT 


BASE.YEAR(I) + J - 1 


return 
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routine CALC.ALT.CODE 
'' Assigns alternative codes to each arc. 
for K = 1 to NUM.NODES (1) 
ALT.CODE(K) = 0 
for I = ( NUM.NODES(1) + 2 ) to N.NODE 
for K = EP(I) to ( EP(I+1) - 1 ) 
ALT.CODE(K) = ( mod.f ( K - EP(I) , N-.ALT ) + 1 ) 


return 


end 


WOON NU &WNY 


routine CALC.ARC 


'' Determines the tail node, life-cycle cost & effectiveness 
'* associated with each arc. 


'' Initialising LCE & Limiting MOE arrays 
for K = 1 to N.ARC 


do 
LCE(K) = 0 
LIM.MOE(K) = 0 
loop 


'' Computations for the present system ( system 1 ). 
‘' All arcs here have the alternative code 0. 
for K = 1 to NUM.NODES (1) 


do 
TAIL(K) = 1 
M = TRANS.YEAR(K+1) - CURR.YEAR 
Lcc(K) = 0.DC * exp.f( O.DC.EF * M ) - O.SV * exp.f( O.SV.EF * M ) 
for I = 1 to NUM.MOE 
do 
MOE.VALUE = O.MOE(I) '' Lcc is made very high 
if ( MOE.VALUE < MIN.MOE.REQ(I,M+1) ) ‘'* iff any MOE require- 
Lcc(K) = 9999999 '' ment level is not met 
LIM.MOE(K) = I 
always 
'' Note : Utility fn. is Value / Max.Value 
if (I= 1 
LCE(K) = LCE(K) + MOE.WT(I) * MOE.VALUE / 180 
else 
LCE(K) = LCE(K) + MOE.WT(I) * MOE.VALUE / 100 
always 
loop 
for J = 1towmM 
do 
AMC = O.MC.P1l * O.MC.P2 * ( J ** ( O.MC.P2 - 1 ) ) 
AOC = 0.0C * exp.f( 0O.OC.EF * J ) 
LCC(K) = LCC(K) + ( AMC + AOC ) * CUM.DF(J) 
for I = 1 to NUM.MOE ‘ 
do 
MOE.VALUE = O.MOE(I) * exp.f( O.MOE.DF(I) * J 
if ( MOE.VALUE < MIN.MOE.REQ(I,J+1) ) '' Checking if MOE 
LCC(K) = 9999999 '* req. were met. 
LIM.MOE(K) = I 
always 
‘*' Note : Utility fn. is Value / Max.Value 
if (I =1) 
LCE(K) = LCE(K) + MOE.WT(I) * MOE.VALUE / 180 
else 
LCE(K) = LCE(K) + MOE.WT(I) * MOE.VALUE / 100 
always 
loop 
loop 
loop 
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‘* Computations for later systems, represented by 
'' arcs with alternative code higher than 0. 


for S = 2 to N. 


SYS 


for J = 3 to N.NODE 
with SYS.CODE(J) = S 


for I = 2 to N.NODE 
with SYS.CODE(I) =S - 1 
do 
M = TRANS.YEAR(J) - TRANS.YEAR(I) - '' Number of years used. 
if ( ( M >= MIN.LIFE(S) ) and ( M <= MAX.LIFE(S) ) ) 
for N = 1 to N.ALT 
do 
TAIL(K) = I '' Start node of the arc K. 
Dey = DC(N,I) * 
exp.f( DC.EF(N,I) * M ) '' Disposal cost & 
SVY = SV(N,I) * '* salvage value in the 
exp.f( SV.EF(N,I) * M) ‘'' last year. 
YR = TRANS.YEAR(I) - CURR.YEAR '' No. of years 
'' from the present. 
LCC(K) = AC(N,I) * CUM.DF(YR+1) + 
( DCY-SVY ) * CUM.DF(YR+1) 
for E = 1 to NUM.MOE 
do 


MOE.VALUE = MOE(N,E, YR+1) 
'' Note : Utility fn. is Value / Max.Value 
if iE = 70) 
LCE(K) = LCE(K) + MOE.WT(E) * MOE.VALUE / 180 


else 
LCE(K) = LCE(K) + MOE.WT(E) * MOE.VALUE / 100 
always 
if ( MOE.VALUE < MIN.MOE.REQ(E,YR#+1) ) 
Lcc(K) = 9999999 '' Checking if MOE 
LIM.MOE(K) = E ‘' req. were met. 
always 
loop 
for Y = 1 toM 


do 
AMC = MC.P1(N,I) * MC.P2(N,I) * '' Annual maint. 
( ¥ ** ( MC.P2(N,1I) - 1) } '' & operating 
AOC = OC(N,I) * exp.f( OC.EF(N,I) * ¥ ) '' costs. 
LCC(K) = LcCC(K) + ( AMC + AOC) * 
CUM. DF (YR+M) 
for E = 1 to NUM.MOE 
do 
MOE.VALUE = MOE(N,E,YR+t1) * 
exp.f ( MOE.DF(N,E,YR+1) * Y ) 
''t Note : Utility fn. is Value / Max.Value 


if ( E=1 

LCE(K) = LCE(K) + MOE.WT(E) * MOE.VALUE / 180 
else 

LCE(K) = LCE(K) + MOE.WT(E) * MOE.VALUE / 100 
always 


if ( MOE.VALUE < MIN.MOE.REQ(E,YR+tY+1) ) 
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111 
112 
113 
114 
115 
116 
ele, 
118 
119 
120 
121 
22 
123 
124 
125 
126 
Pe? 
128 
129 end 


LCC(K) = 9999999 
LIM.MOE(K) = E 
always 
loop 
loop 
compute MAX.LCE as the 
maximum of LCE(K) 
if ( Lec(K) < 9999999 ) 
compute MIN.LCC as the 
minimum of LCC(kK) 
always 
K = K + 1 


loop 


'' Checking if MOE 
'' req. were met. 


MAX.LCE & MIN.LCC are 
needed for computing 

the upper bound on the 
Lagrangian multiplier. 


1 routine ECHO.INPUT 

Pe 

3 '' Echo of important inputs 

4 

5 skip 1 line 

6 print 1 line thus 

7 Sys.Code Min.Life Max.Life Base.Yr. No.Nodes No.Arcs 

8 For I = 1 to N.SYS 

9 print 1 line with I, MIN.LIFE(I), MAX.LIFE(I), BASE.YEAR(I), 
10 NUM.NODES(I), NUM.ARCS(I) thus 

a k** kkk kkk kkk kak KKK 

12 

13 skip 1 line 

14 print 1 line thus 

15 Alternative Name Alt.Code Lead Time 

16 For I = 1 to N.ALT 

Vy print 1 line with ALT.NAME(I), I, LEAD.TIME(I) thus 

18 RAHAKKKKKKKKKKKK KKK kkk kK 

19 

20 skip 1 line 

21 print 2 line thus 

22 Measure of Effectiveness Present System 

23 Code Name Weightage Level Exp. RateofChange 
24 for I = 1 to NUM.MOE 

25 print 1 line with I, MOE.NAME(I), MOE.WT(I), ©O.MOE(I), O.MOE.DF(I) 
26 kkk RAHAEKKKKKKKKKKRE SK wk kkk RAKKKK KKK kaeaek kkk 
27 

28 return 

29 

30 end 
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routine LAGRANGIAN 


'' Solution by Lagrangian relaxation method 


The starting lower bound for the Lagrangian multiplier lamda is 0. 
The upper bound is determined such that it is barely sufficient to 
ensure that the corresponding slope will be positive. Variables 
may be overwhelmed, even when double precision is used, if lamda 
or the cost values are too high. 


0 
( N.NODE * MAX.LCE / MIN.LCC ) + 1 


LAMDA1 
LAMDA2 


skip 1 line 
print 3 line thus 
Iterations of the Lagrangian Relaxation Method 


Lamdal Lamda3 Lamda2 Slopel Slope3 Slope2 Ll L3 L2 


Skip 1 line 
'' Iterations by intercept of supporting planes method. 
‘Iterate' 


LAMDA = LAMDA1 

call BEST. PATH 

call SUM.LCC.LCE 

SUM.LCC.1 = SUM.LCC 

SUM.LCE.1 = SUM.LCE 

SLOPE1 = -1 * ( SUM.LCC.1 - BUDGET ) 

Ll = SUM.LCE.1 - LAMDA1 * ( SUM.LCC.1 - BUDGET ) 


if ( SLOPE1 > 0O ) 
skip 1 line 
print 2 line thus 
The budget is sufficient for cost to be no longer a constraint, 
je. the solution is that of maximum effectiveness over time. 
return 
always 


LAMDA = LAMDA2 

call BEST. PATH 

call SUM.LCC.LCE 

SUM.LCC.2 = SUM.LCC 

SUM.LCE.2 = SUM.LCE 

SLOPE2 = -1 * ( SUM.LCC.2 - BUDGET ) 

L2 = SUM.LCE.2 - LAMDA2 * ( SUM.LCC.2 - BUDGET ) 


if (Slopes. 00) 
skip 1 line 
print 4 line with BUDGET thus 
No solution is possible because of one of the following reasons: 
(1) the given budget of *#kkkAAKA AA JQ Insufficient 
(2) one or more MOE requirement level is too high 
(3) cost values are too high and needs to be scaled down. 
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56 stop 


a7 always 

58 

59 if ( ( SLOPE1 = SLOPE2 ) and ( Lil = L2 ) ) 

60 skip 1 line 

61 print 1 line thus 

62 The cost values are too high for computation & should be scaled down. 
63 stop 

64 always 

65 

66 '' Intercept of the line segments at lamdal & lamda2 yields lamda3. 
67 LAMDA3 = ( SUM.LCE.1 - SUM.LCE.2 ) / ( SUM.LCC.1 - SUM.LCC.2 ) 

68 

69 LAMDA = LAMDA3 

70 call BEST. PATH 

zeal Gall, SUM. LCC. LCE 

72 SLOPE3 = -1 * ( SUM.LCC =- BUDGET ) 

73 L3 = SUM.LCE - LAMDA3 * ( SUM.LCC - BUDGET ) 

74 

Pi) print 1 line with LAMDA1, LAMDA3, LAMDA2, SLOPE1, SLOPE3, SLOPE2, 
76 L1, L3, L2 thus 


T7 KKK aaa KKK KKK KaKRK hhh KKK Kh RK kk hhh kk kkk kh kh ko kK Ok 


79 '' The minima of the piece-wise linear Lagrangian function L 

80 ‘*' is reached, ie. the solution is found, when the intercept 

81 ‘* lamda3 is equal to any 28 the two bounds 

82 if ( LAMDA3 = LAMDA2 ) '' Minima is found & slope there is positive 
83 return ‘' ie. minima itself is the optimal solution 
84 always 

85 if ( LAMDA3 = LAMDAl1 ) ‘t Minima is found but slope is negative. 

86 LAMDA = LAMDA2 '' Lamda2 lies on the solution line segment, 
87 call BEST. PATH 't which is the one with positive slope next 
8B call SUM.LCC.LCE '' to the minima. 

B89 SLOPE3 = -1 * ( SUM.LCC - BUDGET ) 

90 L3 = SUM.LCE - LAMDA3 * ( SUM.LCC - BUDGET ) 

91 print 1 line with LAMDA1, LAMDA3, LAMDA2, SLOPE1, SLOPE3, SLOPE2, 
92 fie ta, £2 thus 


93 kak chk Rak aK khhh kik KKK Ka KKK Kk kkkhkh kk kkk ck kkaK kh kkk & 
94 return 


95 always 

96 

ay '* Lamda3 replaces the redundant bound for the next iteration. 
98 if ( SLOPE3 >= 0 

99 let LAMDA2 = LAMDA3 
100 else 

101 let LAMDA1 = LAMDA3 
102 always 

103 

104 go to ‘Iterate' 

105 

106 end 
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routine BEST. PATH 


'' Forward dynamic programming with criteria function dependent on 
'' the given approach. Provisions are made to ensure that the solu- 
'' tion path has sufficient lead time for each alternative chosen. 


CRIT. FUNC(1) = 0 
for I = 2 to N.NODE '' Initialising crit. func. with a very 
CRIT.FUNC(I) = -99999999 '' low value as O is too high a default 
‘' value to prevent adverse selections. 


for I = 2 to N.NODE 


do 
MAX.CF = -999999999 
for J = EP(I) to (EP(I+1)-1) 
do 


'' Lagrangian relaxation approach to obtain the 
‘' max. eff. for the given budget. 
if ( APPROACH = "Lagrangian" ) 
if ( LCCc(J) < 9999999 ) 
L(J) = ( LCE(J) - LAMDA * Lcc(J) ) 
else 
L(J) = -9999999 
always 
always 


‘' The approach of maximising the cost-eff. 
'' ratio while meeting given requirements for 
'' effectiveness. The budget is ignored. 


if ( APPROACH = "Cost-Eff.Ratio" ) 
if ( Lec(J) < 9999999 ) 
L(J) = LCE(J) / ( LeC(J) - MIN.LCC + 1 ) 
else 
L(J) = -9999999 
always 
always 


'' The approach of satisfying effectiveness 
'* requirements at minimum cost. 
'' yields the lower bound on budget. 


if ( APPROACH = "Satisficing" ) 
L(J) = -1 * Lec(J) 
always 


‘' The approach of maximising effectiveness without 
‘' concern for budget. Yields the upper bound on budget. 
if ( APPROACH = "Max.Eff." ) 
if ( LCC(J) < 9999999 ) 
L(J) = LCE(J) 
else 
L(J) = -9999999 
always 
always 


'' preventing the selection of an arc that cannot 
‘' be used because of insufficient lead time. 


PZ 
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109 
110 end 


ft ALT.CODE(J) > 0 } 
if ( BEST.PRED.NODE(TAIL(J)) = 0 ) 


L(J) 
else 
K = 


if ( K < LEAD.TIME(ALT.CODE(J)) ) 


— ooo 


TRANS. YEAR(TAIL(J)) - 


TRANS. YEAR(BEST. PRED.NODE(TAIL(J) )) 


L(J) = -9999999 
always 
always 
always 
if ( CRIT.FUNC(TAIL(J)) > -9999999 ) 


if ( CRIT.FUNC(TAIL(J)) + L(J) ) > MAX.CF 


MAX. 


always 
always 


loop 
loop 


'' The terminal node with the highest value in the 
'' criteria function yields the solution path that 


CF = CRIT.FUNC(TAIL(J)) + L(J) 
CRIT.FUNC(I) = MAX.CF 
BEST.PRED.NODE(I) = TAIL(J) 
BEST.PRED.ARC(I) = J 


'' is optimal for the given approach. 


MAX.CF = -999999999 
MAX.MAC = ~999999999 
MAX.AAC = -999999999 


for I = 1 to N.NODE 
with SYS.CODE(I) = SYS.CODE(N.NODE) 


x* 


( 
( 


He La) 
INT.DF = 


do 
if ( APPROACH = "Satisficing" ) 
N = TRANS.YEAR(I) ~ CURR.YEAR + 1 
INT.DF = ( ( INT.RATE / 100 ) + 1 ) 
AAC.CF = ( INT.DF * (INT.RATE/100) ) / 
if ( ( CRIT.FUNC(I) / N ) > MAX.MAC ) 
MAX.MAC = CRIT.FUNC(I) / N 
ARG.MAC = I 
always 
if ( ( CRIT.FUNC(I) * AAC.CF ) > MAX.AAC ) 
MAX.AAC = CRIT.FUNC(I) * AAC.CF 
ARG.AAC = I 
always 
always 
if ( CRIT.FUNC(I) >= MAX.CF ) 
MAX.CF = CRIT.FUNC(T) 
ARGMAX.CF = [I 
always 
loop 
return 
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1 ) 


1 routine SUM. LCC.LCE 

2 

3 '' Totals up the life-cycle cost & effectiveness 
4 '' for the solution path. 

5 

6 I = ARGMAX.CF 

7 SUM.LCC = 0 

8 SUM.LCE = 0 

9 until I = 1 

10 do 

fea if ( BEST.PRED.ARC(I) = 0 ) 
12 print 2 lines thus 


13 All paths violates either lead time or MOE level requirements 
14 or both. No feasible solution can be found. 


rs stop 

16 always 

17 SUM.LCC = SUM.LCC + LCC(BEST.PRED.ARC(I) ) 
18 SUM.LCE = SUM.LCE + LCE(BEST.PRED.ARC(I) ) 
19 I = BEST.PRED.NODE (I) 

20 loop 

21 

22 return 

25 


14 


ODNAU AUN 


routine RESULTS. PRINT 
't print-Out of Results 
skip 1 line 


print 6 line thus 
Network of Possible Solutions in Nodal Form 


Note - The criterion function of a node will be very low if all paths 
leading to it fail to satisfy lead time req. and/or minimum MOE levels. 


No. Sys. Year EP Crit.Fn. Best.Pred.Node Best.Pred.Arc Best.Alt. 

skip 1 line 

print 1 line with SYS.CODE(1), TRANS.YEAR(1), EP(1) and 
CRIT.FUNC(1) thus 

a kar kk * hakkk kik 

for I = 2 to N.NODE 
print 1 line with I, SYS.CODE(I), TRANS.YEAR(I), EP(I), 

CRIT.FUNC(I), BEST.PRED.NODE(1), BEST.PRED.ARC(I) and 
ALT.CODE(max.f£(1,BEST.PRED.ARC(I))) thus 
kkk kk kk * kanak kak kik ake kkk 

if ( CRIT.FUNC(ARGMAX.CF) <= -99999 ) 

print 1 line thus 
Warning : Lead-time requirements cannot be met. 

always 

if ( SUM.LCC >= 9999999 ) 
print 1 line thus 

Warning : Minimum MOE level requirements cannot be met. 

always 

if ( APPROACH = “"Satisficing" ) '' There are two possible solutions 
skip 1 line '' for the satisficing approach. 


print 2 line thus 
Minimum Average Annual Cost Solution Path 


( Time value of money is ignored, ie. interest rate is assumed to be 0 ) 


skip 1 line 
print 2 line thus 
Transition Alternative Life-Cycle, Life-Cycle 
Node Year System Employed Cost Effectiveness 
skip 1 line 
I = ARG.MAC 
until I = 1 
do 
print 1 line with I, TRANS.YEAR(I), SYS.CODE(I), 
ALT.CODE(BEST.PRED.ARC(1)), LCC(BEST.PRED.ARC(I)), 
LCE (BEST. PRED.ARC(I)), CRIT.FUNC(I) thus 
zkk* ax * kk KRakKKK kkk RaAKKRK KKK 
I = BEST. PRED. NODE(I) 
loop 
skip 1 line 
print 2 line with INT.RATE thus 
Minimum Annualised Cost Solution Path 


( Monetary interest assumed constant at ***,*** %& per annum ) 


Criterion 
Func.Value 


khkkk kk 
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oa 
94 
95 
96 
oF 
28 
og 
100 
101 
102 
103 
104 
105 
106 
107 


Nod 


ak 


Sol 


Nod 


ax 


skip 1 line 
print 2 line thus 
Transition Alternative Life-Cycle Life-Cycle Criterion 
e Year System Employed Cost Effectiveness Func.Value 
skip 1 line 
I = ARG.AAC 
until I = 1 
do 
print 1 line with I, TRANS.YEAR(I), SYS.CODE(I), 
ALT.CODE (BEST. PRED.ARC(I)), LCC(BEST.PRED.ARC(I)), 
LCE (BEST. PRED.ARC(I)), CRIT.FUNC(I) thus 
x ax 


* ax xzax Kak kak RAKE KA RaaRAK Ka 
I = BEST. PRED.NODE(I) 
loop 
return 
always 


skip 1 line 
print 4 line thus 


ution Path 
Transition Alternative Life-Cycle Life-Cycle Criterion 
e Year System Employed Cost Effectiveness Func.Value 


skip 1 line 
I = ARGMAX.CF 
until I = 1 


do 
print 1 line with I, TRANS.YEAR(I), SYS.CODE(I), 
ALT.CODE(BEST.PRED.ARC(I)), LCC(BEST.PRED.ARC(I)), 
LCE (BEST. PRED.ARC(I)), CRIT.FUNC(I) thus 
x ax x zanx araaaAK Rik RaaKK kak aaah kar 
I = BEST. PRED.NODE(I) 
loop 
skip 1 line 


if ( ( CRIT.FUNC(ARGMAX.CF) > -99999 ) and ( SUM.LCC < 9999999 ) ) 
print 2 line with SUM.LCC and SUM.LCE thus 


The total life-cycle cost and effectiveness for the optimal solution 


are 


ret 


end 


AkaRKKAA AH aNd *htkdkk RHA respectively. 
always 


skip 1 line 

if ( APPROACH = "Lagrangian" ) 
print 1 line with LAMDA3 thus 

Optimal value of the Lagrangian multiplier is *#*kx,arrae 
skip 1 line 

always 


urn 
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APPENDIX B: SAMPLE INPUT FILES 


Stationary. Environment 
3 23 1989 0 


Off 
O-5 
5 10 
5 10 
Alt.A 0 
Alt.B 0 
Weietel td ktddyidii 
ete tek tt ld bia 3 
Lagrangian 
1200 

20°1.3 7 «05 0 0 20 =.2 
100 30 1.3 7 .05 0 0 50 -.2 
POOss 0+ 1.3 7 .05 0 0 50 -.2 
160-530 1.3 7 .05 0 0 50 =-.2 
EOO 30°).3 7 .05)0 O850 -g@2 
100 30 1.3 7 .05 0 O 50 -.2 
HOO 30° 1.3 7 .05 0 0 508-.2 
100 30 1.3 7 .05 0 0 50 -.2 
100 30 1.3 7 .05 0 0 50 ~.2 
HWOge3071.3. 7.05 0 :°0550 =.2 
100 30 1.3 7 .05 0 0 50 -.2 
fodgG30° 1.3.7 .05 0 0 SG) -.2 
LOO 30 1.3 7 .05 6 0 50 -—.2 
100° 30 1.3 7 .05 0 0 Site-.2 
100)30 1.3 7 .05 0 0 50 -.2 
bOoGs30° 1.3 7 .05 0 0 Sde-.2 
rOOesO 1.3 7 .05 0 0 SO =-.2 
Loos3sG 1.3 7 .05 0 0 5G -.2 
JOO 30 1.3 7 .05 0 O 50 -.2 
TOG 30 1.3 7 .05 0 0 50 ~.2 
106°30 1.3 7 .05 0 0 50 =.2 
100230 1.3 7 «05 0 0 507-2 
100°30 1.3 7 .05 0 0 50 =. 2 
10030 1.3 7 .05 0 0 507-22 
HOO 30 1.3 7 .05 0 0 50 -.2 
100 30 1.3 7 .05 0 0 50 -.2 
POO330 1.3.7 .05 0 0 SO =.2 
300 15 1.3 4 .05 0 0 90 -.2° 
300 15 1.3 4' ..05 0 0 90 =—.2 
200) 1590.3 46.05 0 0 $0 =. 2 
300 15 72.3 4).05 0 0790 =-.2 
300 15 1.3.4*°.205 O O0@90 =-.2 
300 15 1.3 4.05 0 0 90 =.2 
200 15 1.3 4 .05 00 90 =.2 
300° °15 1.3 4 .05 0 0 90 =.2 
300 15 1.3 4: .05m@0 0 90 -.2 
300 15 1.3 4 .05 0 0 9ae -E2 
300 15 1.3 §# .05 0 0 SEI -.2 
300 15 1.3,%%05 0 0 3@ -22 
300 15 1.3 4 .05 0 0 90 -.2 
306 15 1.34.05 0 6 30 -.2 
300 15 1.3 @ «05 0°50 90 =.2 
300 15 1.3.45.05 O80 90 —.2 
300 15 1.349.055 O70 90 -.2 
300 15 1.3 41 .05-0,60 90 >= 72 


300 
300 
300 
300 
300 
300 
300 
300 


a 


® 


WWW W WW Ww Ww 
bh bbb hh hb 


Range (nm) 
Accuracy (%) 
Availability(%) 0.3 


95 
25 


80 
80 


70 
70 


120 
120 
120 
120 
120 


140 
140 
140 
140 
140 


95 95 
95.95 


80 
80 


80 
80 


70 
70 


70 
70 


-03 
~03 
-03 
~03 


95 
95 


80 
80 


70 
70 


120 
120 
120 
120 
120 


205 ONOn90 
-05 0 0 90 
-05 0 0 90 
-05 070790 
-05 0 0 90 
-05 0 0 90 
705 OCN0'90 
-05 0 0 90 
0.4 104. 
0.3 87. 
93. 
95  95Ro5e9 > 
95 95/9595 
80 80 80 80 
80 80 80 80 
70 70 70 70 
70 70 70 70 
ae O03 120 - 
=. 035120: = 
aes 12045 
=.03 4220 - 
=~ O03 
s02 95 -.02 
702 95 -=.02 
-02 95 -.02 
-O02 95 -.02 
~ 02 
~O03) 95a. 03 
2037 95) = 203 
w03 95 = 703 
303 6250 =203 
~03 
=. OS. 40 "> 
=. 03 140° = 
=.03.140 > - 
aeO03 1407- 
=2.03 
“02°95 = 502 
-02° 95 =-,02 
P02) 99 = 102 
202 95 =7..02 
-02 
MOOD OD a= 03 
WOs 95 =. 03 
SOomd >. =. 03 
2035955 '=. 03 
as, 


03 


-03 
-03 
.03 


-03 
- 03 
~03 
203 


2D 
35 
35 
25 


NNNNNN NHN NY 


-~.03 
-.02 
= OS 


95 
oS 


80 
80 


70 
70 


25595 
95°95 


80 80 
80 80 


70 70 
70 70 


120 
120 
120 
120 


~02 95 
woes 99 
~02 95 
J02,95 


703°95 
=.03 95 
203° 95 
703 95 


140 
140 
140 
140 


=O02 oo 
-02 95 
-02 95 
2-02 95 


~03"95 
203° 95 
=—,O03— 95 
=..09 Goo 


118 


- 03> 220 
-03 120 
03 120 
~-03 120 


er 8 


-..03 


03 
=-<03 


-.03 140 
-.03 140 
-.03 140 
-.03 140 


. 03 
- 03 
. 03 
~03 


= .03 120 
=. 03 7420 
-.03 120 
-.03 120 
95  =.02 
95 =.02 
95 -.02 
ee eee (0 PZ 
93 = .03 
956e..03 
pe NS ees ONE 
95ai= 208 
-.03 140 
-.03 140 
-.03 140 
-.03 140 
35 «=. 02 
95) =.02 
95=—. 02 
oo amos 
95 =.03 
95°=.035 
99° =208 
95 6a 03 


-.03 


=..03 
.03 
- 03 
203 


Inflation.&.Interest.Rates 
3 23 1989 5 

OLft 

0 5 

5 10 


B99 )599 595 .95 .95 95 
foo. 99299 £95 .95 2.95 
Rea OID 957.95 2.95 
Rooms 9o 6.95 .95 .95 .95 
295 0. 99 


Satisficing 
900 

30.0 1.3 7.0 .05 QO 
mOO-0 31.21.3 7.3 .05 0 
PO4-0 32-4 1.3 7.6 .05 O 
HOG. 2 33.7 1.3 7.9 .05 O 
Mee oe SO. 1.3 8.2 .05 0 
mi7.O 36.5 1.3 8.5 .05 0 
Bas? S820 1.3 8.9 .05 O 
ecw 39.9 1.3 9.2 .05 0 
W31.6 41.1 1.3 9.6 .05 0 
MyiG6-9)42.7 1.3 10.0 .05 0 
142.3 44.4 1.3 10.4 .05 0 
148.0 46.2 1.3 10.8 .05 0 
W53-9 48.01.3 11.2 .05 0 
"60.15 50.0 1.3 11.7 205 0 
mee. on 92.0 1.3 12.1 2.05 0 
i7a22754,.0 1.3 12.6 .05 0 
Mambo o.2 2.3 13.1 .05 0 
Pevoeoo.4 1.3 13.6 .05 0 
194.8 60.8 1.3 14.2 .05 0 
wOez.G Gi.2 1.3 14.7 .05 0 
Zoe oo.? 1.3 15.3 .05 0 
2797 1568.4 1.3 16.0 .05 0 
eaecgeet. t 1.3 16.6 2,05 0 
Zogons7o.9 1.3 17.2 .05 0 
Zao. o776.9 1.3 17.9 .05 0 
eoaeIeo0.0 1.3 18.7 .05 O 
Zoic. o7o9.2 1.3 19.4 .05 0 
woweomto.6 1.3 4.2 .05 0 
Sraz,or16.2 1.3 4.3 2.05 0 
924.5 16.9 1.3 4.5 .05 0 
wove ome.) 1.3 84.7 20S O 
9579-0 18.2 1.3 4.9 .05 0 
Beeeon9.0 1.3 95.1 .05 0 
Wwemorr9.7 1.3 5.3 .05 0 
BeanoreO-.5 1.3 5.5 .05 0 
Aiwmoewe). 3 ko 3 5.7 2.05 0 
di2eeomee- ec 1.3 95.9 .05 0 
aa4cp eS.) 2.3 6.2 .05 0 
g616 24-0 1.3 6.4 .05 0 
aiGeteeo.,0 1.3 6.7 .05 0 
499.5 26.0 1.3 6.9 .05 0 
Boe oe27.O0 1.3 7.2 .05 0 


ooooo0oo0ooocoeoeo0oc“c*deodooo0oooqoo0o00o0o 8 


eoo0o0oqoo°0oeoeo0o0dcodces 


20 

50.0 
52.0 
54.1 
56.2 
58.5 
60.8 
63.2 
65.8 
68.4 
71.2 
74.0 
77.0 
80.1 
83.3 
86.6 
90.0 
93.6 
97.4 
101.3 
105.3 
109.6 
3.9 
Pro. 2 
Zon 2 
Peo 
13320 


90.0 
93.6 
STE S| 
101.2 
105.3 
109% 5 
mi. 9 
118.4 
Me. 2 
128.1 
133.2 
138.6 
144.1 
149.9 
755.9 


! 
e e e ¢* oe oe e®# @#  @  @ @ ee e¢ oe @® @  @  @  @  @ 
NNUNNNN NM NNN NMNNNNNNNNNNNNNN YN 
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NNNNMNNNNMYNNYNNNN 
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540.3. 268.1 1:3 97 35> een G2 es. 2 
561.9 29.2 1.37) 7280s Oso moo Ons. 2 
584.4 30.4 1.3 “8.1 Coe eee!) >. se]. 2 
607.7 31.6 1.3 8.4 32505 07 081c2 32-.2 
632.1 32.9 1.3 8.38.05 Cs0 eco. 6Gy—.2 
657.3 34.2 1.3) [9 057010819722 =.2 
683.65 35..5- . 3) 25 F205 Ont e 05.1 = 2 
711.0°37.0 1.3 “S89 2057080213 245.2 
7139.4,538.4. 1.3 1033) OS807 02221 ee. 2 
769°.0"407071.3 1027) .059030 23077 —.2 
799.8 41.56 71 sale I O50 0 239 59R—. 2 
Range (nm) 0.4 104.6 =,03 

Accuracy (%) OFS ows 2 iO? 

Availability(%) 0.3 93.5 -.03 


95,9595. 95 95.95 95595 95 95595 79sn 
95°95 95 95 95 °95°0'5) 95795 95a 75> Jo 


80 80 80 80 80 80 80 80 80 80 80 80 80 
80 80 80 80 80 80 80 80 80 80 80 80 80 


70 70 70 70 70 70 70 70 70 70 70 70 70 
70 70 70 70 70 70 70 70 70 70 70 70 70 


120 -.03 120 =.03 120 —-.03 120 -.035220>—. 03a Ome 
120 -.03 120 —-.03 120 -.03 120 —.03°5220 ={020Z 0g 
120 —-.03 120 —.03 120 —.03, 120 2e0sel Oe ore Oe 
120 -.03 120 —.03 120 —.03 120 Se03RI 208 Oe One 
120—.. 0375120) —303 


95 =.02 95 —.02°95 —202 995-202 oo 02s eee 
95 =.02 95 -.02 95 =.02 95 =. O24 5te 02a 
95 =.02 95 —-.02 95 —.02 9555.02 895 4-202 eae 
95 -.02 95 -.02 95 —.02 95 —-. 02895 2. Ua oe 
9555-02 95e— 02 


95 =.03 95° -—.03 95 -.03 95 —.03 98> (eee oa 
95 =.03 95 =.03 95 -.03 95 =.03795 72 aGeeoeae ws 
95 =—.03 95 =.03 95 =-.03 95 —.03 95320 n ee oe 
95 =.03 95° -.03 95 —.03 95 —208 09> ene eee 
95 -..03 955.03 


140 -.03 140 -.03 140 -.03 140 -.03 140 —-.03 140 -.03 
140 -.03 140 —.03 140 —-.03 140 —203 1405-70382 OF .0 
140 -.03 140 —-.03 140 —-.03 140 —.035140) 2 o3 Ce. oe 
140 -.03 140 —.03 140 —-.03 140 —2035 14032023408. Oo? 
140 —.03 2460) —-.C¢3 


95 -.02 95 -.02 95 -.02 95 —=.02 507 O 2 eee 
95 -.02 95 -.02 95 -—.02 95 —.02 S52 2G2ee> a 
95 -.02 95 -.02 95 =.02 95 —.02 9554.0 8 oe 
95 -.02 95 -.02 95 3.02 95 —.02 95) 7.07 oe 
95 =—.02 95 so702 


95 -.03 95 =30379S =. 03 99 203 oe Oe oe 
95 =.03 95 =-.03m95 —.03 95 —- 3038383 > ou 
95 -.03 95 =—.03795 —.03 95 —.0355> 6. Uso oa 
95 -.03 95 —.03795)—-.03 95 =303 (956 oe eee 


ey) 


Technological. Change 


See 9 1969 5 


Off 
0 5 
5 10 
5,10 
Alt.A 


A1lt.B 


see he} 
95 
~95 
~95 
ele, 


~95 
~95 
295 
oS 
<99 


0 
0 


oD 
yk, 
~995 
~99 


Satisficing 


650 


100.0 
102.0 
104.4 
106.1 
108.2 
110.4 
112.6 
114.9 
117.2 
19). 5 
i21..9 
124.3 
126.8 
129.4 
ono 
134.6 
£37.3 
140.0 
142.8 
145.7 
148.6 
151.6 
154.6 
157.7 
160.8 
164.1 


300.0 
303.0 
306.0 
309.1 
oie .2 
So. 2 
318.5 
321.6 
324.9 
Bie Br) 
331.4 
334.7 
Boe. 0 
341.4 
344.8 
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799 
295 
=o 
295 


wovuouovooaonnry ws 
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SPANIUONAOWANAKRP HAKAN DWHOADANVUOUNUO AW © 


2 


NWN E&P NUON OW ODO WY 


~95 
oS 
eva) 
oS 


ooooooocooooqe°oeo0od0cdcdcdc0dc@qco000o0 0 000 © 


oooo0o0oo0oocooo0o0o°0o°0o0c dc 0° 


Oo oo oqo oqoooqoqcoqoqooqocqeqcooqoqc°oqcqo coc 3} 


ooooo0oo0ooo0oo0ooo0o09$o;o 


20 

50.0 
51.0 
92.2 
53.0 
54.1 
557.2 
ot ee 
57.4 
58.6 
597 
60.9 
62.1 
63.4 
64.7 
65.9 
67.3 
68.6 
70.0 
71.4 
72.8 
74.3 
75.8 
77.3 
78.8 
80.4 
82.0 


90.0 
90.9 
91". 8 
92.7 
93.7 
94.6 
99.9 
96.5 
97.5 
98.4 
99.4 
100.4 
101.4 
102.4 
Oss 


2 


* 


NONNNNNNNNNNNNNNNNNNNNNNN NN NY 
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NNN MM 1 MM MY YN —| AY 


J 


Kho 
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348.3 28.1 1.3 
S54 6529.2 1.3 
355.3 30.4 1.3 
358.8 31.6 1.3 
362.4 32.9 1.3 
366-1 34.2 1:3 
3695) 35.541 7535 
37364 37.0 1.3 
3722) 38.4 1.3 
380.9 40.0 1.3 
384.7 41.6 1.3 
Range (nm) 
Accuracy (%) 
Availability (%) 
95 95 95 95 95 
95.95 95 95595 
80 80 80 80 80 
80 80 80 80 80 
70 70 70 70 70 
70 70 70 70 70 
1205-03123 7— 
1505-03154 5— 
172 =.03 175 = 
180 —..03 180 — 
180 -=.03 180 .— 
95.07—,02 95.1 
95.67—.02°9 95.7 
96.2 —=.02 36.5 
96.7 —.02 96.7 
96.9 -.02 97.0 
95.07 =..03 95.1 
95.6: "= 7.103: 9 5:27, 
9O.28--035 96.4 
S6./.—.03.96.7 
96.9 =.03 97.0 
140 -.03 144 - 
158 -.03 160 - 
168 =—-03 169 = 
173 =.03 174 .= 
1765 =. 03217 7- — 
95: .0°-2.0220 5.2 
95.6°—.02595 77 
96.2 -.02 96.3 
96.7 —.02,96.7 
96.9 =.02°97-0 
95.0 °=203 995.7 
95.6 —-303 95-7 
96.2 —-,o2 96.2 
96.7 =—.03 96-7 
96.9 —=.03 97.0 


fomd fmd fmt 


80 
80 


70 
70 


.03 


03 
03 


03 


e 


-03 


03 


-03 


03 


-03 


Kmeoownrwnwvowvooodrs 


-02 


-02 


03 


-02 
-02 
-02 
-02 


-03 


.03 


e e e 


PN WO Ore Oh eH OU 
oOo oo ¢ oo © 2 © 
OOaInnoannn vo 


205 


104. 
87. 
93. 


55 


80 
80 


70 
70 


128 
158 
178 
180 


95.2 
95.8 
96. 
96. 


02 
02 


02 


03 


03 
03 


147 
162 


| Neer a 


252 
95.8 
96.4 
96.8 
02 


95-2 
95.8 
96.4 
96.8 


03 
03 


03 


OFOR1 04 oF —. 2 

0 Oel05-5° ].2 

0 0 106.6 -.2 

0 0 107.7 -.2 

0 0 108.7 -.2 

0 0 109.8 -.2 

0-0 13059 -—.2 

O50 sii2.0 =. 2 

0-0-3134) =.2 

O° 0 11453 -.2 

00 115.4 -.2 

6 =. 03 

2 =.02 

5 =,03 

95 95 95 95 95 

95_95-95 735205 

8G 80 80 80 80 
80 80 80 80 80 
70 70 70 70 70 
70 70 70 70 70 

-03 134 -.03 140 

-03 162 -.03 166 

~O3° 179 =. 038760 

-03 180 -.03 180 
=,02°795.3 32702 
=.02 95.9 52.02 
=.02 96.4,—.02 
=.02 96.8 .—.02 
=.03 95.3) —.02 
=.03 9573) =.05 
-.03 96.4 = .03 
-.03 96.8 =2.03 
03.150 =—203" 153 
03 164 -.03 165 
03 171 =.03 472 
03 175 —. 03 yc 
-.02 95.3 >] 902 
=.02 95.9 = 702 
=.02 96.4 "=202 
=.,02 96.68 {2 GzZ 
-,03 95.3 —.02 
=.03 95.9 =.02 
-.03 96.4 -.03 
=-.03 96.87 — 7020 


fz2 


=. 03 
=, 03 
=.03 
-.03 


95.4 
96.0 
96.5 
96.9 


95.4 
96.0 
96.5 
96.9 


=.,03 
—.03 
=. O35 
03 


95.4 
96.0 
96.5 
96.9 


95.4 
96.0 
96.5 
96.9 


145 
169 
180 
180 


=. 03 
=.03 
-.03 
=O 


302 795m 
~O2 96.1 
-02 96.6 
-O2 96.9 


-03 95.5 
~03796.1 
=-03796.c 
=O 9 One 


-.03 
=~.03 
=. 03 
=.03 


i302 
167 
173 
176 


= 02 95. 
-.02 96. 
=202 796. 
-.02 96. 


OAr Ul 


=,03) 957.5 
=.05 96.1 
-.03 96.6 


=.00 96.9 


02 
~02 
-0O2 
nO 2 


03 
-03 
-03 
510) 3) 


02 
202 
- 02 
702 


103 
03 
eens 
= 2.08 


Requirement.Escalation 
3 23 1989 5 

Off 

OF-5 

5 10 

SLO 

Alt.A 0O 

Alt.B 0 


eee 695 295 .95 .95 
moomoo. FID 695 2.95 «695 
oom 5s 95 295 .95 .95 
Boog. 95 7.95 .95 .95 


295° .95 
Lagrangian 
700 

sO-0) 1.3 7.0 5..05 O80 20 = ee 
PoomewsO.3 12-3 7.3 -05 0 @ 50°20 -.2 
meee 3s0.6 1.3 7.6 .05 0 @ 51.0 -.2 
104.4 30.9 1.3 7.9 057000 52225 —.2 
muGet 31.2 1.3 6.2 .05 0 0 53 ae se 
MOGece 31.5 1.3 8.5 .05 0 0 54. -.2 
P1024 31.8 1.3 8.9 .05 0 0 5526 =. 2 
Ieeeoms2.2°1.3 9.2 .05 0 0 56.3 -.2 
Pees ses > 1.3 9.6 .05 0 0 57.4 -.2 
ive esge.6 1.3 10.0 .05 0 0 586.56 -.2 
ieee. tk 1.3 10.4 .05 0 0 59.7 -.2 
ieee s. 5 1.3 10.8 .05 0 0 60.9 -.2 
Weare 3.8 1.3 «11s2 ) .05 0 0 62.1 =.2 
i2eeees4 1 12.3 11.7 .05 0 0 63.4 =—.2 
Pe2o74e54-.9 1.3 12.1 .05 0 0 64.7 -.2 
Moles 34.8 1.3 12.6 .05 0 0 65.9 -.2 
imoeeuso.20l.3 13.1 .05.0 0 67.3 —.2 
Myo s >. odo 139-6 .05 0 0 68.6 =—.2 
MoO o> -9mt. > 14.2 .05 00 70.0 -.2 
miiz2.o 36,20..3 04.7 .05 0 0 71.4 -.2 
M7) 36.6) 8.3 15.3. .05 0 0 72.8 -.2 
mice s7.0 1.3 16-0 .05 00 74.3 -.2 
ieee 37.3 1.3 16.6 .05 0 O TS 6S. 2 
Womoug7. 701.39 17.2 <05 O20 77.3 =.2 
ieee? oe. tl. 3 1749 .05 OF0OF 78.8 = g2 
ioeroeso.ou!.3 L827 .05 0708 80.4 -.2 
Mot 36.9 2.3 1954 .05 OFOF 2.0 -.2 
Ooo! 5.6 1.3 4.2.05 0 0 90.0 =.2 
303.0 16.2 1.3 4.3 O5 (0705 JOR =. 2 
MOoeOmo.9 1.3 4.5 .05 0 0 98 —.2 
aU.) 1.3 42.7 .05 0 0. 9257 —.2 
aize2 1652 1.3 4.9 .05 0 0 9O93%7 —.2 
five 3 1920 1.3 5.1 .05 0 0 S940 6a-.2 
reco 19.7 1.3 5.3 .05 0 0 95257 -.2 
Oemeomzo.5 1.3 5.5 .05 0-0 96.25 -—.2 
S2aeoeet. ot 1. Ss. 5.7.05 0 0) 9725 -.2 
S20.le22.2 1.3 Sao. 2.05 0 O08 9854 -.2 
Soe 4A e259.) 2.3 Gee .05 0 OU O8R4 —.2 
H34.7 24°50 1.3 6.4 .05 0 0100.4 -.2 
aeeOmeaa0 1.3 6.7 .05 0 0 104 =.2 
4245 26.0 1.3 6.9 .05 0 0 102.45 -.2 
Mier omeenOwl. 3 7.2 .05 O 0-103. 54—.2 


123 


348.3 28.12 1.3° 7.5 35057 0n 0g 04-27. 2 
351.8 29.2 1.3 97. 80>) Open 05. 555-2 
355.3 30.4 1.3508. 1965 020R106-67-.2 
358.8 31.16 1.3998.4 4205, 0 20—4 07 7a 2 
362.4 32.9 1.3 8.8 305 OF G08 74-2 
366.1 34.2 1:39. 13305 Osseo ee 2 
369.7 35.5 1.35,925 205 02 0sitOr oe. .2 
373.4 372.0 1435) 929 3205. 030m 2 os. 2 
377.1 38.4 1.3710.3 (205 OSeFii si]. 2 
380.9 40,0 1.3°10.7 20> 0 O73 114735-.2 
384.7 7417.6 1.351101 . 05) OOwerS 74-2 
Range (nm) 0.4 104.6 =.03 
Accuracy (3%) O.3° 5887.2 —FiO2 
Availability(%) 0.3 93.5 -.03 

95 95 95 100 100 100 110 110 110 110 
110 110 110 110 110 110 110 110 110 110 


80 80 80 80 80 
80 80 80 80 80 
70 70 70 72 72 
102 DOD 7 DO 15 
1207 =. 039123 — 
150°=.03 254 <= 
L7i2-. 035175) — 
L805]. Oo eis0y— 
[LeO |. 03 1240)>— 
95.0 =. 029959 
95.00 =.02 95-7 
Sibe2 =.02 96.3 
26.7 —.02 96.7 
965.9. '=.102 (9770 
95.0) =. Ose O) 
hale ar 0/5). SNS U/ 
96.2 = 708996 .3 
96.7 =. Cem 6.7 
96.9 =2303497..0 
140 -.03 144 - 
158) >. 038160 5= 
168 —.03 7169 >— 
173) =. 034A 
1/6 =.03 417 72 — 
95.0 =.025 9551 
95.6 =, 02595 7 
96.2 —.02 95.4 
9647 —-.02)596"7 
96.99 —. 025970 
9520) —20de95cu 
99.6 —-203 795.7 
96.2 —.03 796 
96.7 -.03396—7, 
96.9 —-703897 0 


80 80 80 80 80 80 80 80 
80 80 80 80 80 80 80 80 
72°72. 75°75 75 7 See ao 
75° 75° 7575.75 7537557 > 
03 128 =.03 134 — 7038, 1406 
O3,158 =.05 1262 =~05 7166 
034178 =203. 179 [23052 tse 
O03 180 =.0391260 —-03 160 
03 

=.02°95.2 -.02 953-07 
=-.02 95.8 =2702 OSes «07 
=,02 96.4 —=-.02° 96-4 —.02 
-.02 96.8 -.02 96.83-.07 
=202 

=-.03 95.2 -.03 9535505 
=.03 95.8 —. 03° 95 33-0 
-.03 96.4 —.03 9674 — 03 
-.03 96.8 =.03 96735-7702 
ar ONE. 

O03 147 =.03 150 =—2037 152 
O3 162 =.03 164 —.034 6s 
03.170 =.03 471 =.033ve2 
035175 =.03 175) -- 0s 
03 

=.02 95.2 =-.02 9s nee. 
=.02 95.8 =-.,02 95° oe 
-.02 96.4 -.02 96,4 7243C@2 
-.02 96.8 -.02 96.852 7072 
=.02 

=.-03 95.2 —.03 95.358 —104 
=.03 95.8 =. 03 257582 04 
=.03 96.4 -.03 9G747_—.02 
=.03 96.8 —.03 962707 —.08 
=0'5 


24 


210 1105120 
110° 1768110 


=. 03 
= 203 
= 03 
= 30s 


95.4 
96.0 
96.5 
S659 


95.4 
96.0 
96% 5 
26n 


= 10s 
=.03 
=.03 
=. 03 


95.4 
96.0 
96.5 
96.9 


95.4 
Sern (0, 
oer S 


96.9 


=03 
=O 
=.03 
See OS 


145 
169 
180 
180 


=702"°95 30 
—BO2 96.1 
=302 96.6 
oc 9629 


= 3 “95.5 
=—s03 96.1 
=.03 96.6 
por 01S) oN ce, 


=e O18) 
= 2035 
= 203 
=, 03 


love, 
167 
173 
176 


—m02 95.5 
woe) 96. 1 
aa 25 9 © 316 
a2 96. 9 


ee 0)5) Behe 
-03 96. 
-O3 96. 
=. 03 96. 


Onr WN 


=O 


=. 03 


=~ O02 
nF 


=.103 
=.03 


Upgrading.Scheduling 
3 23 1989 5 

Off 

0 5 

5 10 
5210 
Alt.A 0O 

Alt.B 0 

-95 .95 .95 .95 .95 
95 .95 .95 .95 .95 
mo wom .95 .95 .95 
wo> .95> .95 .95 .95 
s95 .95 


Cost-Eff.Ratio 


780.08 
20-0 1.3 7.0 
moog; 0OF30.3 1.3 7.3 
82.0 30.6 1.3 7.6 
¥O04.4 30.9 1.3 7.9 
MOG.r S1-2 1.3 8.2 
woe2 91.5 1.3 8.5 
110.4 31.8 1.3 8.9 
M2.6°32.2 1.3 9.2 
Pome 2255 1:3 9.6 
117.2 32.8 1.3 10.0 
P99. > 33.1 1.3 10.4 
S329) 55.5 1.3 10.8 
Wass 393.8 1.3 11.2 
ieio.6 34.1 1.3 11.7 
WOle4 34.5 1.3 12.1 
ieieoeoao 1.3 12.6 
anos oo.2 1.3 13.1 
Og soo.) 1.3 13.6 
140.0 35.9 1.3 14.2 
awe 3Oe2 14.3 14.7 
ies! 26.0 1.3°15.3 
148.6 37.0 1.3 16.0 
fo 397.39 1.3 16.6 
Peieeeo 97.7 1.3.17.2 
toe S62.) 1.3 17.9 
Hoe 36.5 1.3 18.7 
128.1 38.9 1.3 19.4 
900.0 15.6 1.3 4.2 
93-0 16.2 1.3 4.3 
906.0 16.9 1.3 4.5 
woo.) 17-5 1.3 4.7 
Me. 2 28-2 1.3 4.9 
ao. 3 19,0 12.3 5.1 
Seo. toe7 1.3 5.3 
Seoegoo 1.3 5.5 
waaeo 21.3 1.3 5.7 
320,72) 22.2 1.3 5.9 
aot 4 236) 1.3 62 
334.7 24.0 1.3 G4 
m8 ,0, 25.0 1.3 ° 6.7 
341.4 26.0 1.3 6.9 
ga4-8527.0 1.3 7.2 


GN 
es: 


-95 


o-oo OO Oo oO 0 0 Oooo Oo Oo 2d Ooo Oo oo © 
Ooqoooo0ooqoqocceoocooqoqoqcoqcooqoqoqcocqdo } 


ooeoo0ooo0ocoooo0oo0oo0o0coe 


Oooqocqooqooqoqoqoqo°qco cd 


Of Owodnw Dreher OnwWOnN PrP OAIN A WNFON O 


Nf&b ll LuyuyranrwsnWdwo 


NNNNNNNNNNNNNNNNNNNNNNNNN N VY 


* 


NNNNNNNNNNNNN NN 


— 
bo 


Cr 


3483520. 1215 
35256) 29 22 01 
3553 30.4003 
358.8531. 6.1-2 
362.4. 32. 901.3 
366. 1934.2 0173 
369 27 8395). 54.0 
373.4 37.20 41.3 
377 21636 24 
380.9 40,051.35 
384.7 4166-13 
Range (nm) 
Accuracy (%) 
Availability (%) 
95 95 95 100 
110 110 110 110 
80 80 80 80 80 
80 80 80 80 80 
70°70 70. 72572 
1D “2S: 7547557) 
i205. 032s 
150 3=.03>1545— 
L723. 03 oe 
180 -.03 180 = 
180 =, 033 lS05— 
95.0 -. 02-95. 
SS..6 =s02 S57 
SiGe 2 = 02 es Oe 
96.7 =.02796.7 
Bagel SO 97). 10, 
95.0 =20g5 95. i 
95.6°—,03595.7 
96.2 -=. OS 876, 3 
96.7 —-.03 96.7 
96.9 —703"97-.20 
140 203 91445— 
158 =262 1607 — 
166 =2-G2) 169 —— 
173 -=208. 174 <— 
176 -.03 177 - 
95.0 =—2027 95.1 
95.6 =.02 95.7 
96.2 "—702eg6.3 
96.7 =.02) 96-7 
96.9 —.02597.0 
95.0 =,03 295-8 
95.6 =. 09mg or 
96.2 —.03596.3 
9657 =. 03 ea6n7 
96.97 —, 0399720 


7.52057 0"05 101755 —. 2 

71 <8 2OSmORO Ot oo 2 

8.1 2.05 -O305106.65——-2 

8.4 ..05 OCs0e 07.7 ae 

&.8 05707 OFlOety ye —-2 

9.1 OS2080 VCoe es. 2 

9.59.05. USO. Aor <= ..2 

9.9.05 7070-11220 =.2 
102.37 3095070 2 ial 2 
10.7 2.05050 a3 see 
Li.1 305700 (1 ioe4 ae 

0.4 104.6 =2038 

0.3 88722 = oe 

O23 935 =203 

100 100 1707 110-1160 srie 

110 110 110 110 110 110 
80 80 80 80 80 80 80 80 
80 80 80 80 80 80 80 80 
V2 72 73 15°79 7587 cao 
75 75 75 75275) 7 5 oe 
-03 128 -.03 134 -.03 140 
03-158 =.03 162 =.63, 166 
OS e178 — 7035 a7 9” — 7031 so 
Os 180 -—.03 180 —303 “180 
03 

=.02°95.2 -. 02795. 56-7 
=.02 9558 —.02) 99. 3a 
=-.02 96.4 =.027 S6e44-.07 
=,.02 96.8 —.02596.0 2-07 
= 02 

=-.03 95.2 —503495. 28s 
=.03 95.8 —-2.038795 79s. 
-.03 96.4 -.03 96.452 .63 
=-03 96.8 =-.03 96.ey oc 
=, 03 

O3 147 -—.03 150 —203yIs3 
O3 162 —-.03 164 =.03e%6> 
O03 170 =-.03 17) -—JOsmHG 2 
OS 175 =-.03 175 =. 035 va 
73 

=.02° 95.2 =502995 33a 
=-,02 95.8. —.02- 95793 — ur 
=.02 96.4 =—.02)96.45—202 
=.02° 96.8 =.02 9G) ene 
=—02 

= 0395 2-2 (= .0 30 9 Sesser 
=.03595.28 —-.03995.2—4,—.0.5 
=.039796.4 —.03°96.45— 0c 
=O3" 96.5 —,.05—596 78.05 
= 03 


126 


110 110 110 
170 BWiGe ire 


—-.03 
= 203 
= 0s 
OS 


95.4 
96.0 
9625 
96-9 


9554 
96.0 
96.5 
96.9 


=.02 
=. 08 
=O 
a5) 6) 


95.4 
960 
965 
96.9 


95.4 
26.0 
96.5 
96.3 


145 
169 
180 
180 


sear OS 
=. O35 
=. Os 
—.08 


=O) 6 Ou 
=O2796..1 
=.'027.96..6 
=U 27.9 62.0 


= OBES o> 
tO S eo Oren! 
=.03' 916.6 
=.03 596.9 


195 
£67 
73 
176 


cal O 
=.03 
= OS 
eo | 


=——O2 Goi. D 
Oe oo. 1 
— 20287 6.6 
=202 4596.9 


—- Ome o . 9 
5 (3)2) Silene 
=—~ 0396.6 
=—.05 96.9 


= 20 
=.02 
=202 
8b 


=O2 
O02 
Oe 
;02 


703 
.03 
-03 
203 


Development. Lead.Time 
Bee. os 1969 5 

On 

0 5 

5 10 

5 10 

Alt.A 1 

Alt.B 0O 


Peo 699 «9S .95 95 
moo. foo 295 2.95 .95 
f99.99 .95 .95 .95 .95 
Jomo ooo 1.95" .95 .95 


7-99 6995 
Lagrangian 
634.52 
30.07) .3a97.0 .05 0 0 20 Ss 
MOCO O0. 3 1.36 7.3 .05.0 0 SQR0.—.2 
B2.2.0 30.6 lee 7.6 .05 0 0 SiO -—.2 
Maan sO.9 12.3 7.9 .05 0 0 52.2 =—.2 
meGe lt giv? 1.3 8.2 .05 0 0 5350.—.2 
BG-2 J4mo 1.3 §8.5 .05 070 S54 -.2 
20.4) Sieo 1.3 8.9 .05 ©.0 552.2 —.2 
2 eorce2-2 1.3 “9.2 .05° 0970 5623 — 72 
BOmgues 2. Oo 1s 30 39.6 .O05 8050 57.4e—..2 
ive o2.0 1.3 7°80.0 .05 0 0 58.6 —22 
ies. > JS.) 1.3 10.4 .05 0-0 59277 -—.2 
w3eoes5.o 1.3 10.8 .05 O86 60.9 =.2 
izaeseoseo) 3.3 11.2 .05 Ow 62.1 =—22 
Pzoeees4.t 1.3 11.7 ~.05 0 0 6324 -.2 
Ot. S34. 3 12.1 405 0-0) 647 =a 
elo) 4. 3 12.6 2.05 0 07 65.9 -—F2 
ma Ome. e ol. kd 205 0:10 267.3 —.2 
BO ego. o elo 13s.6 6605: 0 U0! 6 6mG68..UG =. 2 
PaO. 0us5.9 153 14.2 .05 0 0 70.0 -.2 
Pace osO.2c 5 14.7 .05 0.0. 74 =. 2 
is? SO.00te Set. > ~05 00°) 7288 =—.2 
rae.6 37.081,3 716.0 05 OBO) 7485 =.2 
mol.6 37.3 Voo0l6.6 ~.05 O80” 7on8 -—82 
m20nG S727 229 017.2 .05 OO 7753 =< 72 
mie eooes 125 17.9 .05 0 0 7E58 —.2 
MoO0eGsG.5 1.3 18.7 .05 0 0 8054 -.2 
Gate 98.9 V2.9 19.4 ..05 OF0. 82n0 —.2 
S800-0°515.6 1.3 4.2 .05 © 0 9G20 =—.2 
303.0 16.2 1.3 4:3 .<05 0 0 SQ@zo9 =. 2 
BOG .OrrG.9 1-3 4.5 .05 0 0 @2.5 —.2 
309-17. oo) 467 «05 0.0 92,7 =.2 
312.2 1652 1.23,.4.9 2.05 0-0 98.7 =-.2 
go. 19.0 7leog yo. 1 205 0 0 9aG =—22 
Seno 19.7 o3ay 5.3 305 0 OF 95E5 =.2 
eeomeO.)> les 5.5 205 0 O0- 9685 —.2 
eames «=©655l 7 OD 0 0 9725 =.2 
O20 ege.e 1.3 555.9 .05 0 Ommeot.4 -—.2 
D314 25% 1 eo 2 205 O Ome 994 —.2 
334-7 24.0 123% 6.4 205 0 Omp0CE4 —.2 
336.0 25.0 1.3 6.7 2.05 0 0-104 -.2 
341-4 26.0 1.3 6:9 .05 0 0 102-4 -.2 
44627. 0 leon 7.2.05 0° 0 20825 -—.2 


P20) 


34872357 28-15 1-2 
351.8 29.2 1.3 
355.3 30.4 1.3 
358-6 31.681 .3 
362.4 32.90)..3 
366.1 34.201 .3 
369.7 35.5981 .3 
373.4 37.0 1.3 
377.1 38.4 1.3 
380.9 40.0 1.3 
384.7 41.6 1.3 
Range (nm) 
Accuracy (%) 
Availability (%) 
95° 9595-100 
110 110 110 110 
8O 80 80 80 80 
80 80 80 80 80 
70° 70 70 J2072 
75 75 75 75575 
1207 —..035 123 5— 
150 —.03°> 154, -. 
72 =-203 175... 
180 -.03 180 -. 
180 —.03 18C —; 
95.0 —.02 95.1 
95.6 =.02 095.7 
96.2 =.02 96.3 
96.7 —-.02 96.7 
96.9 5— 202,977.06 
95.0 —.03 95.1 
95265-2037 99.7 
96.2 -.03 96.3 
965.7 9=,.0 39 396.07 
$6.9 —,03 9720 
140 -.03 144 -. 
158 —-.0@3 160 —- 
168 —.03 169) — 
173 °=2.037 174 2. 
176 —.03 177 —. 
95.07-—.02,95-1 
95.67 —=2027 95777 
96.2 —.02796.3 
96:27 =—202 96.7 
96.9 -.02 97.0 
95.0 .-203> 95.1 
95.6 —.03° 95.7 
96.2 =. 057 96-35 
96.7 =3039 96.7 
96.9" =—303,°9 7.20 


110 110 110 
110 110 110 


7.5 ~05 0 0 104-5 
7.8 205 0 0 105.5 
8.1 .05 0 0 106.6 
‘8.4 .05 0 O 107.7 
8.8 .05 0 O 108.7 
9.1 205-0 071058 
9.5 205 7070s 0-9 
9.9 2.05,°0 O7112.0 
10.32.05 07 0-1ts 1 
10.7 2.05 0 (08114 3 
11.1 .05 0 O 115.4 
0.4 104.6 =. 09 
0.3 87.2 —.0Z 
0.3 993. =.03 
100 100 110 
110 110 110 
80 80 80 80 80 80 
80 80 80 80 80 80 
72 72 75°19 15.75 
78 75 75 75 75 75 
03 128 —.03 134 -—. 
O39 158>—2035 162 —-: 
03 178 —-.03 5179 = 
O32. 180 -—.0375180 = 
03 
=.02 95.2 5=.02 952 
=.02 95.8 =.02 95= 
-.02 96.4 —-.02 96. 
—-.02 96.8 -.02 96. 
-.02 
=203, 95.2 —=203 952 
=.03 95.8 —.0327 95 
-~.03 96.4 -.03 96. 
|=.03 96.8 —.03> 96. 
=.03 
O3n 1475 —..03 
~03 162 —.03 1643—. 
Os. 170 —.03 17) 
03° 275 =—.035175 ae 
03 
=.02 95.2 [202 95e 
=-.02 95.8 =.02 952 
=.02 96.4 —.,02 9G. 
-.02 96.8 —.02 96; 
=~ 02 
=.03), 95.2 =. 0305. 
=.03°95.8 —.103)°95. 
=.03 96.4 =—.03 9Ge 
=.03 96.8 —.03° 9Ge 
=.03 


80 80 
80 80 


7575 
75575 


03 140 
03 166 
03 180 
03 180 


3 =202 
9 =.02 
4 -.02 
8 -.02 


3 =203 
9 =.02 
4 -.03 
8 —03 


1507=.03 253 


03 165 
03 172 
03 176 
3 =. 02 
9 =. 02 
4 -.02 
8 -.02 
3 =.03 
9 =. 08 
4 =.03 
8 =-.09 


128 


110 110 110 
110 110 110 


- 03 
o03 
=.03 
- 03 


95.4 
96.0 
IGF 2 
96.9 


95.4 
36.0 
96.5 
96.9 


-03 
=.09 
=. 03 
=-. 02 


95.4 
96.0 
56.10 
96.9 


95.4 
96.0 
96.5 
96.9 


145 =<.03 
G69 203 
180 -.03 
T80 -.03 


-702 95.5 
=-.02 96.1 
-.02 96.6 
-.02 96.9 


=-.03 95.5 
=203° 96.1 
=-2703 96.6 
-.03 96.9 


=.03 
=—03 
=03 
= 20 3 


155 
167 
173 
176 


02 95.9 
=-02 96.1 
-.02 96.6 
= 202 96.9 


—-.03 95.5 
-.0@3 96.1 
=.3 96.6 
=702 96.9 


~O02 
=e02 
=.02 
O72 


03 
-.03 
03 
03 


Sau2 
-.02 
-.02 
-.02 


-.03 
— 03 
-.03 
=—s03 


APPENDIX C: SAMPLE OUTPUT FILES 


Dynamic End-Of-Life-Cycle Planning Model 
Case of Technological Change 
Approach - Lagrangian Budget = §$ 601.74 


Note : Budget is used only for the Lagrangian relaxation approach and is 
ignored in the other approaches. 


Iterations of the Lagrangian Relaxation Method 


Lamdal Lamda3 Lamda2 Slopel Slope3 Slope2 Ll L3 L2 
O. ~025 2eUGns) —3CG709 9-32-79 27.04 18.8 11.4 66.7 
wUZ5 -029 2.080 oer -.01 27.04 11.4 11.4 66.7 
-029 -032 2.080 meet) 27.04 27.04 11.4 11.4 66.7 
noes 2032 7032 -.01 27.04 27.04 17.4 dae4 L124 


Network of Possible Solutions in Nodal Form 


Note - The criterion function of a node will be very low if all paths 
leading to it fail to satisfy lead time req. and/or minimum MOE levels. 


No. Sys. Year EP Crit.Fn. Best.Pred.Node Best.Pred.Arc Best.Alt. 


iL 0 1989 0 0. 

2 eos 9 i 1.411 1 1 0 
3 i 1990 2 -649 1 Pa 0 
4d ao oa 3 =". > ial 1 3 0 
5 1932 4 <1493,7-4 1 4 0 
6 171993 5-1.000E+007 1 5 0 
a 1 1994 6-1.000E+007 1 6 0 
8 2 1994 7 -4.107 2 q 1 
2 2 1995 ) -4.476 3 11 1 
10 2 1996 13 =5 012 4 17 1 
11 er oo7 19 = 3.659 5 25 al 
12 2 1998 27 =6..502 5 33 a 
13 201999 37 -6.876 5 43 i} 
14 2 2000 49 =7 oe 5 55 1 
15 222001 59 =. 0OCEs0G7 a 60 2 
16 2 2002 67-1.000E+007 5 67 al 
17 2 2003 73-1.000E+008 0 0 0 
18 2 2004 77-1. 000E7008 0 0 0 
19 3) ANNES, 73 =7 2.772 8 79 1 
20 3 2000 81 =) ie 2 83 1 
21 32001 85 =72.962 10 89 1 
22 a7 2002 pe -8.285 11 97 i] 
Zo 322008 9S =O. og r2 107 iI 
24 3 2004 109 -8.860 13 119 1 
25 a 2005); 127 =9..075 14 ileal 1 
26 S22006. 133 =9.3 16 14 141 1 
ou 322007 145 =-9,.510 14 151 1 
28 Baezoos 157 =97..657 14 161 1 
Zo 3.2009 169 =97 55 14 oly ant 1 
30 3 2010 181-1.000E+007 14 182 2 


31 3 2011 191-1.000E+008 LS 191 1 
32 3 2012 199-1.000E87008 16 IETS) 1 
33 3 2013 205-1.000E1008 0 0 0 
34 3 2014 209-1.000E+608 0 0 0 


Solution Path 


Transition Alternative Life-Cycle Life-Cycle Criterion 
Node Year System Employed Cost Effectiveness Func.Value 
19 1999 3 1 Add 2 OO 5.029 —J,7¢2 
8 1994 2 1 321.559 4.702 -4.107 
2 1989 1 6) -20.000 ay Be 1.411 
The total life-cycle cost and effectiveness for the optimal solution 
are 574.70 and 10.506 respectively. 
Optimal value of the Lagrangian multiplier is .03180 


130 


Dynamic End-Of-Life-Cycle Planning Model 
Case of Upgrading.Scheduling 
Approach - Satisficing Budget = $ 780.08 


Note : Budget is used only for the Lagrangian relaxation approach and is 
ignored in the other approaches. 


Sys.Code Min.Life Max.Life Base.Yr. No.Nodes No.Arcs 


al 0 5 1989 6 6 
2 5 10 1994 11 72 
| 5 10 Wispsie 16 132 
Alternative Name Alt.Code Lead Time 
Alt.A 1 0 
Alt.B Ai 0 
Measure of Effectiveness Present System 
Code Name Weightage Level Exp.RateOfChange 
1 Range(nm) -400 104.600 OO 
2 Accuracy (%) .300 87.200 -.020 
a Availability (3 .300 93.500 -.030 


Candidate Arcs For Solution Path 


Note - LCC is made to be unfavorably high if the minimum requirement 
level for any component MOE of LCE has not been met. The limiting MOE 
is given as the highest index of the MOE(s) that failed to meet the 
minimum requirement level(s) during the life-span covered by the arc. 


Arc Alt. Tail Head Life-Cycle Limiting 

No. Code Node Node Effectiveness Cost MOE 
1 6) 1 2 «1 F45 —-20.00 0 
2 6) 1 3 1.5287 27.67 6) 
3 0 1 4 2<2632 80.95 0 
4 0 1 5 2.97839999999.00 1 
5 O 1 6 3.67489999999.00 1 
6 0 1 7 4.35309999999.00 2 
w 1 2 8 4.7024 302.34 0 
8 2 2 8 4.9501 385.58 0) 
9 1 2 9 5.41609999999.00 1 
10 2 2 9 5.7008 415.63 0 
11 1 3 9 4.74309999999.00 1 
12 2 a 9 5.0034 374.09 0 
13 1 2 10 6.11099999999.00 1 
14 2 2 10 6.4318 443.29 0 
15 1 3 10 5.46269999999.00 al 
16 2 3 10 5.7617 403.70 0 
17 1 4 10 4.8083 299.01 0 
18 2 4 10 5.0436 3162.25 6) 
19 1 2 a 6.78779999999.00 3 
20 2 2 11 7.14369999999.00 3 
2i 1 3 ii 6.16359999999.00 1 
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NRMP NRF NEP NHEPNEP NEP NEP NEPNPNPNEFPNPNPNEP NP NPNPNPNEPNEP NEP NEPNEPNPNPNEPNHEN 


SN ODNNYNNAKHDUMNUNNYNNAAHDONNKHGP& HSI IAAHAONNKHFLAWWINIAARAUUHHWWNNAARAUMUS BWWNN UU & Ww 


6.5004 430.97 
5.53779999999700 


So5083 391.28 
4.8860 270.43 
5.0842 Sod 2 


7.44699999999.00 
7.836799899999" 00 
6.84609999999.00 
7.21969999999'. 00 
6.24819999999.00 


SG. 5529 418.01 
B02 71 313.96 
5.8550 399.75 
4.9637 278.45 
5.1247 340.96 


8.08899999999.00 
8.51169999999700 
7.51089999999200 
7.92009999999. 00 
6. 939999352700. 00 
7.21//99992299 00 
6.34879999999.00 


6.6054 406.06 
5.7164 320.43 
5.3076 369209 
5.0290 269.00 
iG Jhon) 330.40 


8.15829999999-00 
87°60209999999 .00 
7261369999999 206 
JO S3 B99 959999. Oe 
7205159995999 00 
1233029999999). 00 


6.4494 359-376 
6.6580 395-7035 
5.4 309 7) 
5.9340 357299 


8. 26979999999 .00 
&8.67139999999.00 
7.73599999999.00 
8.04779999999.00 


7.1632 396.23 
7.3944 418.86 
6.5340 347.41 
6.6945 383.44 


8.40249999999.00 
8.74059999999.00 
765627999999. 0C6 
8. 159999999 -06 
7.2500 382.67 
7.4349 406.84 
8.53509999999.00 
8..609893999999 00 
7296109999999. 00 
Si pbeys chee ee ieiele 6 0 


132 
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JY 
118 
119 
1 20 
P20 
122 
123 
124 
A Pate, 
126 
127 
128 
122 
130 
138 


PNP NPN EP NUPNPNPNPNENHPNHPNEPNEPNPNHEPNEPNEFPNEPNEPNEPNEPNEPNEPNEPNHPNEP NEP NH 


=a) 
OUMMOOUUMABDWOWUAWADS ~ 


be 
H OO 


11 


B264629999999.00 
$205799999999. 00 


5.0290 
5.1529 
5.7905 
5.9340 
5.0943 
5.1934 
6.5340 
6.6945 
5.8667 
5.9806 
5.1472 
5.2216 
12.2970 
7.4349 
6.6187 
6.7470 
5.9275 
6.0130 
5.2002 
5.2497 


259.72 
337.32 
296.82 
365.87 
263.82 
327.16 
335.28 
392.24 
303.43 
355.21 
254.57 
a1 261 
371.05 
416.48 
340.50 
381.12 
292.74 
345.21 
229.69 
308.53 


7. 96109977999. 00 
8515592799999. 00 


223509 
7.4932 
6.6872 
6.7835 
5.9884 
6.0454 
D2 Oo 
eat, ) 


374.97 
404.96 
328.45 
370.79 
266.54 
Joe 72 
237.36 
299.26 


8.64659999999.00 
8.85799999999 .00 
8.06389999999.00 
Bo219799999997. 00 


7.4269 
i. O207 
6.7558 
6.8200 
6.0493 
6.0778 
5.3060 
5.2936 


361.66 
39422 
701.03 
360.88 
272.86 
J22-99 
229.22 
290.87 


8.75809999999.00 
8. 9273 97999-00 
8.14729999999.00 
8.26419999999.00 


Loose 
7.5742 
6.8244 
6.8565 
6.1102 
6.0959 
5.3466 
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333.510 
384.04 
306.08 
350.70 
263.45 
a7 19 
206.34 
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132 
138 
134 
135 
136 
37 
138 
132 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
uBsyt 
152 
53 
154 
15 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
7 
V2 
173 
174 
175 
Ls 
177 
178 
179 
180 
181 
182 
183 
184 
185 
186 


NRPNEFPNP NEP NEPNPNEFPNEPNEP NH NP NHNEFPNENENPNEPNEPNEPNEPNEPNEP NP NUP NEP NFP NEF N 


Das2lo 282.34 
8.84849999999.00 
GF 7527 999999'.100 
$2230 579999997. 00 
&.30847999999.-00 


7.57392 336.97 
7.6146 373.47 
6.8930 295.48 
6.8769 341.56 
6.1568 239.38 
6.1283 308.21 
5.3692 213.78 
5.3376 274.17 


8 .93889999999.00 
9.02329999999.00 
8.31389999999.00 
8 .35289999999.00 


726558 eee 7 
7.6373 364.02 
6.9455 270.30 
6.9134 332.17 
6.2034 245.62 
6.1464 299.60 
5.4277 206.71 
5. 3543 266.52 


9202929999999 -00 
9.07159999999 .06 
§ 239719999999. 00 
$233.7 7779999999°..00 


7.7134 229.05 
7.6778 354.26 
6.9981 275.42 
6.9339 323.17 
6.2500 237.47 
6.1646 291256 
5.4682 185.59 
>. 209) 258.95 
9.11969999999.00 
9.09849999999.00 
8.46109999999.00 
8.42209999999.00 
tet i 303.13 
7.1005 344.91 
7.0506 266.25 
6.9543 314.78 
6.2967 215.24 
6.1827 283.58 
5.4806 192.99 
5.3815 251277 


9.18889999999 00 
9.14667999999'.00 
8 .524699999929-00 
8.44699999999.00 
7.8299 293.02 
7.igse 2396.21 
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16:7 1 a7, 30 7.1031 242.97 0 
188 2 ney 30 6.9748 306.43 0 
189 1 18 20 6.3109 221.65 0 
190 2 18 30 6.1970 276.03 0 
191 1 BS 31 9..258199599999700 3 
m2 2 15 31 pe Wad Ey ae Bae eee SENS Ci, 3 
193 u 16 31 8.58879999999.00 3 
194 2 16 31 8.47179999999.00 3 
iS '5 1 a7 sial 7.8882 268.76 0 
196 2 17 a 7.7459 327032 0 
97 a 18 ol 7 ebio2 248.45 0 
198 2 18 al 6.9908 298.54 0 
IMS 2, i 16 32 9.32739999999..00 3 
200 2 16 32 9. 20059999999 .00 3 
201 1 17 32 8.65259999999.00 3 
202 2 17 32 8.49669999999.00 3 
203 1 18 32 7.9060 273.39 0 
204 2 18 32 Ia G37 319.33 0 
205 1 17 33 S227 662799999. 00 3 
206 2 17 33 S227 39999999.00 3 
207 AL 18 33 8.67209999999.00 3 
208 2 18 33 8251619999999. 00 3 
209 i 18 34 9.41779999999.00 2 
210 2 18 34 927248609999999.00 3 


Network of Possible Solutions in Nodal Form 


Note - The criterion function of a node will be very low if all paths 
leading to it fail to satisfy lead time req. and/or minimum MOE levels. 


No. Sys. Year EP Crit.Fn. Best.Pred.Node Best.Pred.Arc Best.Alt. 


al 0 1989 0 0. 

2 2 1989 1 20.000 1 1 0 
3 i199 0 2 “27.6006 1 2 0 
4 loo | 3 =00. 950 1 3 0 
5 ees D2 4-1.000E+007 ] 4 0 
6 1 bie he 7} 5-1.000E+007 a 5 0 
a JL TRS BS | 6-1. 00051007 1 6 0 
8 2 1994 7 =282).339 2 7 1 
9 2 1995 S) =395. 6355 2 10 2 
10 21996 13 -379.964 4 a7, 1 
11 2 1997 1s, -458.632 3 22 2 
12 2 1998 27 -498 .963 4 32 2 
3 as ASS ks 37-1.000E+007 2 37 1 
14 2 2000 49-1 .000E+007 3 49 1 
15 2eao0l 59-1 .000E+007 4 og 1 
16 2 2002 67-1.000E+008 0 0 0 
17 2 2003 73-1, 000E7008 0 0 0 
18 2 2004 77-1. 00GE+006 0 0 0 
19 Oo hisyele) ie, =538.O068 8 79 1 
20 3 2000 81 =57 9% 163 8 81 1 
21 3 2001 85 =6172622 8 85 a 
22 a 2002 oF =653.389 8 a1 i 
23 3 2003 99 -708.410 10 103 1 


24 
25 
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27 
28 
Zo 
30 
at 
32 
33 
34 
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2005 
2006 
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24 
10 
4 


Minimum Annualised Cost Solution Path 


Transition 
Year 


2004 
1996 
1991 


109 

i2t 
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145-1. 
15721. 
16952. 
18it=i. 
19 Te. 
199-1. 
205e1. 
209-1. 


Annual Cost Solution Path 
( Time value of money is ignored, 


System 


3 
2 
1 


-741.623 
=791 7333 
-835 -Ye2 
OOOE+007 
OOOE+007 
OOOE+007 
OOOE+008 
OOOE+008 
OOOE+008 
OOOE+008 
OOOE+008 


Alternative 
Employed 


1 
1 
0 


ie. 


E33 
125 
137 
149 
157 
169 


oo0o0 0 
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interest rate is assumed to be 0O ) 


Life-Cycle 
Cost 


3617-659 
299.014 
80.950 


( Monetary interest assumed constant at 


Node 


22 
8 
2 


Transition 


Year 


2002 
1994 
eyes, 


System 


3 
2 
1 


Alternative 
Employed 


1 
1 
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Life-Cycle 
Cost 


371.056 
3022329 
-20.000 


Life-Cycle 
Effectiveness 


7.427 
4.808 
2.263 


§.000 % per annum ) 


Life-Cycle 
Effectiveness 


7.257 
4.702 
~775 


Criterion 
Func. Value 


-741.623 
-379.964 
-80.950 


Criterion 
Func. Value 


—-6537 350 
=282,. 329 
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